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Preface 



My lectures on the effective field theory for heavy quarks, an expansion around the 
static limit, concentrate on the motivation and formulation of HQET, its renormal- 
ization and discretization. This provides the basis for understanding that and how 
this effective theory can be formulated fully non-perturbativcly in the QCD coupling, 
while by the very nature of an effective field theory, it is perturbative in the expan- 
sion parameter 1/m. After the couplings in the effective theory have been determined, 
the result at a certain order in 1/m is unique up to higher order terms in 1/m. In 
particular the continuum limit of the lattice regularized theory exists and leaves no 
trace of how it was regularized. In other words, the theory yields an asymptotic ex- 
pansion of the QCD observables in 1 jm - as usual in a quantum field theory modified 
by powers of logarithms. None of these properties has been shown rigorously (e.g. to 
all orders in perturbation theory) but perturbative computations and recently also 
non-perturbative lattice results give strong support to this "standard wisdom". 

A subtle issue is that a theoretically consistent formulation of the theory is only 
possible through a non-perturbative matching of its parameters with QCD at finite 
values of 1/m (Sect. [4~4] ) . As a consequence one finds immediately that the splitting 
of a result for a certain observable into, for example, lowest order and first order 
is ambiguous. Depending on how the matching between effective theory and QCD is 
done, a first order contribution may vanish and appear instead in the lowest order. For 
example, the often cited phenomenological HQET parameters A and Ai lack a unique 
non-perturbative definition. But this does not affect the precision of the asymptotic 
expansion in 1/m. The final result for an observable is correct up to order (l/m) n+1 
if the theory was treated including (1/m)™ terms. 

Clearly, the weakest point of HQET is that it intrinsically is an expansion. In 
practise, carrying it out non-perturbatively beyond the order 1/m will be very difficult. 
In this context two observations are relevant. First, the expansion parameter for HQET 
applied to B-physics is AQCD/m-b ~ l/(rom.b) = 1/10 and indeed recent computations 
of 1/m.b corrections showed them to be very small. Second, since HQET yields the 
asymptotic expansion of QCD, it becomes more and more accurate the larger the mass 
is. It can therefore be used to constrain the large mass behavior of QCD computations 
done at finite, varying, quark masses. At some point, computers and computational 
strategies will be sufficient to simulate with lattice spacings which are small enough 
for a relativistic b-quark. One would then like to understand the full mass-behavior 
of observables and a combination of HQET and relativistic QCD will again be most 



useful. Already now, there is a strategy (de Divitiis et al., 2003 a de Divitiis et al, 



20036 Guazzini et ai, 2008), which is related to the one discussed in Sect. 5.3 and 



which, in its final version combines HQET and QCD in such a manner. For a short 
review of this aspect I refer to (Tantalo, 2008 1 . 
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Introduction 



1.1 Conventions 

Our conventions for gauge fields, lattice derivatives etc. are summarized in the ap- 
pendix. 

1.2 The role of HQET 

This school focuses on lattice gauge theories. How does heavy quark effective theory 
(HQET) fit into it? The first part of the answer is that HQET is expected to provide 
the true asymptotic expansion of quantities in powers (accompanied by logarithms) 
of 1/m, the mass of the heavy quark, with all other scales held fixed. The accessible 
quantities are energies, matrix elements and Euclidean correlation functions with a 
single heavy (valence) quark, while all other quarks are light. A full understanding of 
QCD should contain this kinematical region. 

The second part of the answer has to do with the challenge we are facing when 
we perform a Monte Carlo (MC) evaluation of the QCD path integral. This becomes 
apparent by considering the scales which are relevant for QCD. For low energy QCD 
and flavor physics excluding the top-quark, they range from 

pa 140 MeV over m D = 2 GeV to m B = 5 GeV. 
In addition, the ultraviolet cutoff of Auv — cT 1 of the discretized theory has to be 
large compared to all physical energy scales if the theory discretized with a lattice 
spacing a is to be an approximation to a continuum. Finally, the linear extent of space 
time has to be restricted to a finite value L in a numerical treatment: there is an 
infrared cutoff Together the following constraints have to be satisfied. 

Air = <c mn , ... ,mu,ms <C a -1 = A uv (1.1) 

The infrared and the ultraviolet eff ects are system atic errors which have to be con- 
trolled. Infrared effects behave as (Liischer, 19861 0(e~ Lni7r ) and are known from 



chiral perturbation theory (Colangelo et al, 2005) to be at the percent level when 
L > 4/m^ pa 6 fin, while the UV, discretization, errors are O((ora quar k) 2 ) in 0(a)- 
improved theories^] With a charm quark mass of around 1 GeV we have a requirement 



of a < l/(2m c ) . . . l/(4m c ) pa 0.1 ... 0.05 fm ( |Kurth and Sommer, 2002[ ) and thus 



L/« « 60 ... 120. (1.2) 

x See Peter Weisz' lectures for the general discussion of discretization errors and improvement of 
lattice gauge theories. 
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Including b-quarks would increase the already rather intimidating estimate of L/a 
by a factor 4. It is thus mandatory to resort to an effective theory where degrees of 
freedom with energy scales around the b-quark mass and higher are summarized in the 
coefficients of terms in the effective Lagrangian. A precise treatment of this theory has 
become very relevant because the search for physics beyond the Standard Model in 
the impressive first generation of B-physics flavor experiments has been unsuccessful 
so far. New physics contributions are very small and even higher precision is needed 
both in experiment and in theory to possibly reveal them. HQET is a very important 
ingredient in this effort. 

Before we focus on our topic let us note that a factor two or so in L/a may be saved 
by working at somewhat higher pion mass and extrapolating with chiral perturbation 
theory, see M. Golterman's lectures. 

1.3 On continuum HQET 

1.3.1 Idea 

We consider hadrons with a single very heavy quark, e.g. a B-meson. Physical intuition 
tells us that these will be similar to a hydrogen atom with the analogy 

hydrogen atom : heavy proton + light electron 

B-meson : heavy b-quark + light anti-quark 

b-baryons : heavy b-quark + two light quarks 



When we take the limit m = mb — > oo ( "static" ) the b-quark is at rest in the 
rest-frame of the b-hadron (B, Ab, . . . ). In this situation, we should be able to find 
an effective Lagrangian describing the dynamics of the light quarks and glue with the 
heavy quark just representing a color source. Corrections in 1/mb should be systemat- 
ically included in a series expansion in that variable. The Lagrangian is then expected 
to be given as a series in Dk/m where the covariant derivatives act on the heavy quark 
field and correspond to its spatial momenta in the rest-frame of the heavy hadron. 

Before proceeding to a heuristic derivation of the effective field theory, let us note 
some general properties of what we are actually seeking, comparing to other familiar ef- 
fective field theories. In contrast to the low energy effective field theory for electroweak 
interactions, where the heavy particles (W- and Z-boson, top quark) are removed com- 
pletely from the Lagrangian we here want to consider processes with b-quarks in initial 
and/or final states. The b-quark field is thus contained in the Lagrangian and we have 
to find its relevant modes to be keptrl 



Another important effective field theory to compare to is the chiral effective theory, 
covered here by Maarten Golterman. Main differences arc that this is a fully relativis- 
ts theory with loops of the (pseudo-) Goldstone bosons and that the interaction of 
the fields in the effective Lagrangian disappears for zero momentum. The theory can 
therefore be evaluated perturbatively. It is also called chiral perturbation theory. In 

2 However, when one carries out the expansion to include 1/m 2 terms, also a whole set of terms 
generated by b-quark loops in QCD which do not contain the b-quark field in the effective theory 
have to be taken into account. An example are 4-fermion operators made of the light quarks, just as 
they appear when one "integrates out" the W and Z-bosons in the Standard Model. 



and so on. 
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contrast, the b-quarks in HQET still interact non-perturbatively with the light quarks 
and gluons. This effective field theory therefore needs a lattice implementation in order 
to come to predictions beyond those that can be read off from its symmetries. 

1.3.2 Derivation of the form of the effective field theory: FTW trafo 

Strategy 

Our strategy is to carry out the following steps, which we discuss in more detail below. 

• We start from a Euclidean action. 

• We identify the dominant degrees of freedom for the kinematical situation we are 
interested in: the "large" components of the b-quark field for the quark and the 
"small" components for the anti-quark. 

• We decouple large components and small components, order by order in Dk/m 

[ tp^D k /m ij\ <C V'hV'h ]• This assumes smooth gauge (and other) fields. It is 
thus essentially a classical derivation. The decoupling is achieved by a sequence of 



Fouldy Wouthuyscn-Tani (FTW) transformations (see e.g. (Itzykson and Zuber, 



1980)), following essentially (Korner and Thompson, 1991). 



• The irrelevant modes are dropped from the theory (often it is said they are in- 
tegrated out) . Their effects are not expected to change the form of the local La- 
grangian, but just to renormalize its parameters. Still it could be that local terms 
allowed by the symmetries happen to vanish in the classical theory. Thus the sym- 
metries have to be considered and all terms of the proper dimension compatible 
with the symmetries have to be taken into account. 

• At tree level the values of the parameters in the effective Lagrangian are given 
by the FTW transformation. In general (i.e. for any value of the QCD coupling) 
they have to be determined by matching to QCD: one expands QCD correlation 
functions in and compares to HQET. This part of the strategy will be 
discussed in detail in later sections. 

Identifying the degrees of freedom 

We consider the free propagator of a Dirac-fermion in Euclidean space, in the time / space- 
momentum representation^] 



S(x ;-k) = J d 3 xe- lkx (^(x)V(0)) = J ^-f k ° x ° [ihtf* 
= S+(x ;k) + S-(x ;k) , 



(1.3) 



with 



S + (x ;p) = d(x )-^-e- E( ^°P + (u), P + (u) = 1 , Ufl = pjm 

MP) 1 



(1.4) 



S_( ao ;p) - fl(-*o)^ye £ (P)-P_(") , P-(u) = 1± f^ 



3 The expectation value (.} refers to the Euclidean path integral, here with the free Dirac action. 
We suggest to verify these formulae as an exercise. 
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where p^ is the on-shell momentum, i.e. 

po = iE(p) = i\/m 2 + p 2 . (1.5) 

Here S+(xq; p) describes the propagation of a quark from time t — to i = xq and 
S-^oj p) describes the propagation of an anti-quark from t = —xq to t = 0. Since the 
Euclidean 4-velocity vector u satisfies u 2 — u^u^ = —1, the matrices P £ {P+,P_} 
are projection operators, 

[P{u)f = P(u) , P + (u)P_(u) = 0, P+(«)+P_(«) = l. (1.6) 

They allow us to project onto the on-shell components of a quark with velocity u. 

The "large" field components corresponding to the quark are given by the projec- 
tion 

^M(aO = P(uMx) , ^ u (x) = 4>(x)P{u) (1.7) 
and the "small" ones, the anti-quark field, are 

lM*0 = P{-u)i>{x) , fo )U (x) - ^(a?)P(-«) , (1.8) 
such that for free quarks 

d 3 xe-*P x (Vvu(^vJ0)) = S+(x ; p) (1.9) 



and similarly for the anti-quark(^] 

For a b-hadron with velocity u, the fields iph.u(x), u (x) are expected to be the 
relevant ones with the other field-components giving subdominant contributions in 
the path integral representation of correlation functions (or scattering amplitudes in 
Minkowski space), while for a b-hadron ip^ u( x )i V'fi u( x ) are expected to dominate. 

In the presence of a gauge field 

When a gauge field is present, we therefore expect an effective Lagrangian for the 
b-hadrons in terms of iph,u, V'h « P lUS a term for the anti-quark. When we rewrite the 
Dirac Lagrangian in terms of these fields, 

£f= Jj(m + V)iP (1.10) 

there are mixed contributions which involve 

V x =i^, = (£ M „ + UftUu) D v , (1.11) 

where the derivative is projected orthogonal to u„. Analogously we have 

2?||=7„Djl, Dl = -u„D v Uu. (1.12) 

From our general consideration of the kinematical situation that we want to describe, 
acting on the heavy quark field is to be considered small (compared to to). In 

4 The terms "large" and "small" components are common ly used when discussing the non- 
relativistic limit of the Dirac equation for bound states, see e.g. |Itzykson and Zuber, 1980| . 
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contrast, £)ji applied to the field will yield approximately = u^m. We therefore 
carry out an expansion with 

Vuip = O(m) V, 

(1.13) 

= 0(1) tp 

and all other fields, such as F^ v , treated as order one. This is often called the power 
counting scheme. 

FTW trafo and Lagrangian at zero velocity 

Having identified the expansion, we perform a field rotation (FTW transformation) 
to decouple large and small components order by order in 1/m. First we consider the 
special case of zero velocity, 



u k = : T>\\= L>o7o , = D k j k , 

The FTW transformation is 

i> -> V' = c s tA , S= ±-D klk = -S+ , 



(1.14) 



(1.15) 



Its Jacobian is one. The Lagrangian written in terms of the transformed fields, 

Sf=i/{V + m)ip' , (1.16) 
yields a Dirac operator (note that S acts to the right everywhere) 

V + m = c- s (V + m)e- s . (1.17) 
Expanding C - s = 1 - 5 + \S 2 - . . . in S = 0(l/m) yields 

P' + m = P + m + {-S,P + m}+±{-S,{-S,2? + m}} + ... (1.18) 

" v ' ^ v ' v v ' 

O(m) o(l) 0(l/m) 

In the evaluation of the different terms we count all fields and derivatives of fields 
(e.g. Fpi,) as order one except for Dq acting onto the heavy quark field. We work out 
the expansion up to order 1/m. A little algebra yields 

V + m + {-S, V + m} = A)7o - ^[7fe7o^feo + -fkl F kl + 2D k D k ] (1.19) 

with o~nv = § hv, 7„] , F fe; = [Dfc, A] and 

i{-5,{-5,C + m}} = - 1 -[ T a H F w + 2 J D feJ D fe ], (1.20) 
2 v ^ 4m i 

-D fc7fc +0(l/m) 



6 Introduction 



such that 

V' = A)7o-^[ IkloFko + ^a kl F kl +D k D k } + 0(l/m 2 ). (1.21) 




In the static part, -Do7o, the large and small components are decoupled, but one of 
the 1/m terms, 7fc7o-FfcO: is off-diagonal with respect to this split. We therefore seek a 
second transformation ip" — e s tp' to cancel also that term, namely we want 

{-S',V' + m} = -L 7fc7o F fco + 0(l/m 2 ) . (1.22) 
2m 

The simple choice S' = j^jolkFko does the job. Now we have the classical HQET 
Lagrangian 

^h stat + + K°* + + °& (1-23) 

^ h stat =^ h (m + A>)<Ah, P + V> h = V<h, ^ h P+=^ h , F± = ^ (1.24) 
^ff at =^ E (m-D )V £ , P_Vh=^E, ^eP-=? E , (1-25) 

^ h (1) = -{Oun + Ospin) , JSff > = -((5 kin + O sp in) , (1-26) 

correct up to terms of order 1/m . We introduced 

O kia (x) = ^ h (a:) D 2 ^ h (x) , O spia (x) = ^ h (x) <r • B(a:) tfh(a:) , (1.27) 
O kin {x) = fc(x) D 2 ^h(^) , O spin (x) = fc(x) cr ■ B(x) fa{x) , (1.28) 
o'k = 2 e ijk°~ij 7 B k = i\eij k Fij , (1-29) 

and the heavy quark fields are the transformed ones, i.e. we renamed tpu — > iph etc. 

Depending on the process/correlation function, just the heavy quark part or just 
the heavy anti-quark part of the Lagrangian will contribute, but there are also pro- 
cesses such as B — B oscillations where both are needed. 

It is worth summarizing some issues that arose in this formal derivation. 

• Assuming D k = 0(1) means that this is a classical derivation: in the quantum 
field theory path integral we integrate over rough fields, i.e. there are arbitrarily 
large derivatives. 

As emphasized before we therefore take this as a classical Lagrangian. Its renor- 
malization will be discussed later, guided by dimensional counting. 

• The derivation is perturbative in 1/m, order by order. This is all that we want. 
In this way we expect to obtain the asymptotic expansion in powers of 1/m. 

• We note that there are alternative ways to derive the form of the Lagrangian. One 
may integrate out the components ip^'&h m a path integral and then perform a 
formal expansion of the resulting non-local action for the remaining fields in terms 



of a series of local operators (Mannel et at, 1992). Another option is to perform a 



hopping parameter expansion of the Wilson-Dirac lattice propagator. The leading 
term gives the propagator of the static action; see exercise |1.1| 
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FTW transformation and Lagrangian at finite velocity 
At finite velocity the transformation is given again by eq. (1.15) but with S 



D^7 M /(2m). For the lowest order (static) approximation, just the anti-commutator 

{V^^} = \{D^Dl}25^ + \[Dl l ,Dl][ ltl , lv ] 

is needed. Since Z)j| = = D^D^ and the second term just involves a commutator 
of derivatives, we see that {T)^,^} = 0(1). Consequently we find 

2= + Z>")^M + ^Ph,u( m + vl )^Ku + 0(l/m) 

= ^h,u{ m - iu n D tt)ih,u + ^h,u( m + iu » D n)' l Ph,u + 0(l/m) (1.30) 
with the projected fields eq. (1.7) and eq. ( |l.8| ). pj 



Let us add a few comments on the finite velocity theory, since we will not discuss it 
further. 

• 0(4) (or Lorentz) invariance is broken. One therefore has to expect a different 



renormalization of Dq and (or as is usually said, a renormalization of u (Chris- 



tensen et al, 20001 Mandula and Ogilvie, 1998)). 



The operator —iD^Uk is unbounded from below. Since it enters the Hamiltonian 
the theory seems to contain states with arbitrarily large negative energies. Result- 
ing problems in the Euclidean formulation of the theory have been discussed in 
the literature ( Aglietti et al, 1992, Aglictti, 1994), but a compelling formulation of 



the theory seems not to have been found. There are also no modern applications 
of the finite velocity theory on the lattice. We will therefore concentrate entirely 
on zero velocity HQET from now on. 

1.3.3 Propagator and Symmetries 

The continuum propagator. 

We consider the static approximation at zero velocity and the latter always from now 
on. The static Dirac operator for the quark is just Do +m so its Green function, Gh, 
(the propagator) in a gauge field A^{x) then satisfies 

(d xa + A (x) + m)G h (x, y) = S(x - y) P+ . (1.31) 

The solution of this equation is simply 

Gh(x,y) = 0(x - 2/ ) exp(-m(x - y ))'Pexp (- / dz A (z , x) 1 S(x - y) P+ , 



i/ii 



(1.32) 



were V denotes path ordering (fields at the end of the integration path to the left) . In 
the same way the propagator for the anti-quark is |^] 



5 One ca n also obtain this Lagrangian by performing a boost of the zero velocity theory I Horgan 
|ei al., 2 009). In the quoted reference also the next to leading order terms are found. 

6Note Pexp {- /*° d 2o A ( 2o ,x)} = Pexp {- J™ dz A (z , x)} f . 
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G h( x ,y) = %o - x ) exp(-m (y - x ))Vexp 



f dz A (z ,x)X5(x~y)P- 
Jya ) 



(-d X0 - A (x) + m)G^{x, y) = 5{x -y)P-. (1.33) 

The mass appears in a trivial way, with an explicit factor exp(— m \x — for any 
gauge field A„. This exponential decay is then present also after path integration over 
the gauge fields in any 2-point function with a heavy quark, 

C h (x, y; to) = C h (x, y; 0) exp(-m (x - y )) . (1-34) 

An explicit example is 

C^ p (x, y; to) = (^(xfrsiphix) ijj h (y)j5My)) > C 1 - 35 ) 



with ip\( x ) a light-quark fermion field. Eq. ( 1.34) means that to shifts all energies in the 



sector of the Hilbert space with a single heavy quark (or anti-quark). We may remove 
to from the effective Lagrangian and add it to the energies later. We only have to be 



careful that to > in cq. (1.31 1, cq. (1.33) selects the forward/backward propagation. 
Therefore we set 

^=^ h (D Q + e)i> h , J?i™=^(-D + e)^, = Effi + m , (1.36) 

where the limit e — > + is to be understood. 

We note that after performing this shift of the energies, there is no difference in 
the Lagrangian of a charm or a b-quark if both are treated at the lowest order in this 
expansion. We turn to discussing this as well as other symmetries of the static theory. 

Symmetries 
1. Flavor 

If there are F heavy quarks, we just add a corresponding flavor index and use a 
notation 

i>h -> V^h = (V>hi,-..,V'hF) T > V'h -> V> h = bPhv ■ ■ >^hir) ( L37 ) 
if h stat = MD + e)Vh ■ (1.38) 

Then we obviously have the symmetry 

M^-tVfaix), Mx)^Mx)V\ VeSU(F) (1.39) 

and the same for the anti-quarks. Note that this symmetry emerges in the large mass 
limit irrespective of how the limit is taken. For example we may take (F = 2 with the 
first heavy flavor identified with charm and the second with beauty) 

TOb — to c = c x Aqcd j or TOb/TO c = c , TOb — > oo (1-40) 



with either c or d fixed when taking TOb — > oo. 
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2. Spin 

We further note that for each field there are also the two spin components but the 
Lagrangian contains no spin-dependent intera ction. The associated SU(2) rotations 
are generated by the spin matrices eq. (1.291 (remember that i/>h , "0h are kept as 
4-component fields with 2 components vanishing) 



1 



o k 

CTfc 



(1.41) 



where the symbol a k is used at the same time for the Pauli matrices and the 4x4 
matrix. We here are in the Dirac representation where 



7o 



1 

-1 



1 





1 



(1.42) 



The spin rotation is then 

Mx) -> e iQfcff * Vh(x) , Mx) ^Mx)e~ la ^\ (1.43) 
with arbitrary real parameters It acts on each flavor component of the field. Ob- 



viously, the symmetry is even bigger. We can take V € SU(2F) in eq. (1.39). This 
plays a role in heavy meson ChPT (Wise, 1992 Grinstein et at, 1992 Burdman and 



Donoghue, 1992). 



3. Local Flavor-number 

The static Lagrangian contains no space derivative. The transformation 
tfh(a:) e l " (x) Vh(z) , M*) -> Mx)e~^ , 



(1.44) 



is therefore a symmetry for any local phase r)(n). For every point x there is a corre- 
sponding Noether charge 

Qh(x) = i(j h (x)ip h (x) [= %p h (x)~f Q ip h (x)} (1.45) 

which we call local quark number. It is conserved, 

doQ h (x) =0Vx. (1.46) 

1.3.4 Renormalizability of the static theory 

Our effective field theory is in the category of local field theories with a Lagrangian 
made up from local fields. In d space-time dimensions, standard wisdom says that such 
theories are renormalizable if the mass-dimension of the fields in the Lagrangian does 
not exceed d. Ultraviolet divergences can then be absorbed by adding a complete set of 
(composite) local fields with mass dimension smaller or equal to d to the Lagrangian. 

According to this (unproverj^]) rule, the static theory is renormalizable. The possible 
counter-terms have to share the symmetries of the bare Lagrangian. They are easily 



7 Power counting as discussed by Peter Weisz at this school is not applicable here, since the prop- 
agator does not fall off with all momentum components. 
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found. From the kinetic term in the Lagrangian eq. (1.361 we see that the dimension 
of the fields is [V>h] =3/2. Only 2-fermion terms with up to one derivative are then 



possible. Space-derivatives are excluded by the local phase invariance eq. (1.44). We 
then have the total quantum Lagrangian 



JzfhO) = CiOi(iz) + c 2 2 (x) 

Ox(x) = ip h (x)iph.(x) , 2 (x) = ip h (x)D ip h (x) , 



(1.47) 
(1.48) 



where the convention c 2 = 1 can be chosen since it only fixes the unphysical field 
normalization, and c\ = Sm has mass dimension [Sm] = 1 and corresponds to an 
additive mass renormalization. From dimensional analysis and neglecting for simplicity 
the masses of the light quarks, it can be written as Sm — (eig$ + e 2 g^ + . . .) A cut in 
terms of the bare gauge coupling go and a cutoff A cut , which in lattice regularization is 
A cu t = l/ a - For a static quark there is of course no chiral symmetry to forbid additive 
mass renormalization. 

This is the complete static Lagrangian. After the standard QCD renormalization of 
coupling and light quark masses, all divergences can be absorbed in Sm, i.e. an energy 
shift. Flavor symmetry tells us that with several heavy flavors, Sm is proportional to 
the unit matrix in flavor space. Energies of any state are then 



jtiQCD j^stat 



Sm - 



Sm— 



m = El)f 



Sm=0 



(1.49) 



Here mbare and Sm compensate the linear divergence (self energy) of the static theory, 
while m is finite. Note that there is no symmetry which would suggest a natural way 
of splitting mbarc into Sm and m. This split is arbitrary and convention dependent. 
The quantity Sm is often called the residual mass. 

A rigorous proof of renormalizability to all orders in perturbation theory has not 
been given but we note the following. 

• Perturbative computations have confirmed the standard wisdom. These computa- 



tions reach up to three loops in dimensional regularization ( Chetyrkin and Grozin, 



2003 Grozin et al, 2008), while in various different lattice regularizations 1-loop 



computations have been carried out 


Eichten and Hill, 1990a Eichten and Hill, 


1990c Eichten and Hill, 1990S Boucaud et al., 19891 Boucaud et al, 1993 


Flynn 


et al, 1991 ( Borrclli and Pittori, 1992[ 


Kurth and Sommer, 2001 ) Kurth anc 


Som- 


mer, 2002, Delia Morte et al, 2005|Pal 


ombi, 2008,Guazzini et al, 2007||Grimbach 


et al, 2008HPalombi et al, 2006[|Blossier et al, 2006|) 



• We will see non-perturbative results which again yield a rather strong confirma- 
tion. 

• Nevertheless a proof of renormalizability would be very desirable. 

1.3.5 Normalization of states, scaling of decay constants 

For the discussion of the mass-dependence of matrix elements we have to think about 
the normalization of states. Standard, relativistic invariant, normalization of bosonic 
one-particle states is 



(p\p') Iol = (2n) 3 2E(p)S(p-p') 



(1.50) 
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The states have a mass-dimension [ |p) rc i ] = — 1. The factor E(p) introduces a spurious 
mass-dependence. In the large mass limit, relativistic invariance plays no role and we 
should choose a mass-independent normalization instead. The standard convention for 
such a non-relativistic normalization is 

(p|p')N R = <p|p') = 2 (2tt) 3 *(p - p') (1.51) 

with [|p)] = -3/2 and 

|p)rci = VE(P)\P)- (1-52) 



Consider as an example where the normalization of states plays a role, the leptonic 
decay of a B-meson, B~ — > t~v t . The transition amplitude A for this decay is given 
to a good approximation in terms of the effective weak Hamiltonian. It factorizes into 
a leptonic and a hadronic part as 

A cx (t u\t(x)^(1 - j 5 )D T (x)\0) (0|u(x)7^(l - j 5 )b(x)\B-) . (1.53) 
Using parity and Lorentz invariance, the hadronic part is 

(0\u(x)^(l-y 5 )b(x)\B-(p)) = (0\A^x)\B-(p)) = Pt ,f B e ipx (1.54) 
in terms of the flavored axial current 

= u{x)~i^-y b b{x) . (1.55) 

There is a single hadronic parameter /b (matrix element) parameterizing the bound 
state dynamics in this decay. We note that it is very relevant for the phenomenological 
analysis of the CKM matrix ( Antonelli et al, 2009[ ) . 

We may now use HQET to find the asymptotic mass-dependence of /b for large 
m = nib- Since to lowest order in \jm the FTW transformation is trivial, the HQET 
current is just 

A* QET (x) = < at (x) + 0(l/m) , < at (x) = u(xhol 5 ^) ■ (1-56) 

The static current Ag tat has no explicit mass dependence. In static approximation we 
then have 

(0|v4 s o tat (0)|B-(p = 0)) = $ stat , (1.57) 
with a mass- independent $ stat . Its relation to /b, 

$ stat = m B 1/2 Po / B = m^fB , (1.58) 



takes eq. (1.521 into account (j>q — E(0) — me). We arrive at the prediction 

/b = + 0(l/m b ) , ^ = 3^1 + 0(l/m c ) . (1.59) 

V m B JD V m B 



The latter use of course assumes Aqcd/toc *C 1. We will see later that these predictions 
are modified by the rcnormalization of the effective theory. 



12 Introduction 



1.3.6 HQET and phenomenology 

Heavy quark spin/flavor symmetry is very useful to classify the spectrum in terms of 
a few non-perturbative parameters or predict relations between different masses, e.g. 

mg, — m B 1=3 m D , — m|, , (1.60) 

771 B' — 771b ~ 771 rj' — 777D , (1-61) 

where 777,3* , tod* are the vector meson masses and with tob' , ttid' we indicate the first 
excitation in the pseudo-scalar sector. The first of these relations has been seen to be 
approximately realized in nature. 

More detailed statements about semi-leptonic transitions B — > Dlv, B* D*lv 
are possible. In the heavy quark limit for both beauty and charm these are described 



by a single form factor, the Isgur Wise function, instead of several (Isgur and Wise, 



1989 Isgur and Wise, 


1990k 


reviews, e.g. ( 


Neubert, 


19941 



helps to understand lattice results for states with a b-quark. 



Exercise 1.1 Static quarks from the hopping parameter expansion 

Consider a Wilson quark propagator in a gauge background field. Evaluate the leading non- 
vanishing term in the hopping parameter expansion (with non-zero time-separation). Check 
that it is the continuum HQET propagator (restricted to the lattice points) up to an energy 
shift. Even if this is a nice piece of confirmation, note that one here takes the limit k — > 
corresponding to ma — > oo, while the true limit for relating QCD observables $^ CD to those 
of HQET is 

$ HQET ~ lim lim$ QCD , 

m—t-oc a — >0 

in that order! 



2 

Lattice formulation 



We start with the static approximation. The 1/m terms will be added after a discussion 
of the renormalization of the static theory. 



2.1 Lattice action 



For a static quark there is no chiral symmetry. Since we want to avoid doublers, we 
discretize a la Wilson (with r = 1). The continuum D iph(x) is transcribed to the 
lattice as 



A)7o H( V o + Voho - aVSVo} , 
and with P+i/'h = V'h ; P-^Ph = V'Rj we have the lattice identities 
D iph(x) = Voiph(x) , D Q ^5(2) = VoV'h^) • 



(2.1) 



(2.2) 



For later convenience we insert a specific normalization factor, defining the static 
lattice Lagrangians 



Jzfh 



1 



1 + aSm 
1 

1 + aSm 



^> h (x)[Vo + 5m]ip h (x) , 



(2.3) 
(2.4) 



The following points are worth noting. 

• Formally, this is just a one-dimensional Wilson fermion replicated for all space 
points x, see also exercise 

• As a consequence there are no doubler modes. 

• The construction of a positive hermitian transfer matrix for Wilson fermions 
(Luscher, 1977 Montvay and Minister, 1994) can just be taken over. 

• The choice of the backward derivative for the quark and the forward derivative 
for the anti-quark is selected by the Wilson term. We will see that this selects 
forward/ backward propagation and an e-prescription as in eq. (1.36 ) is not needed. 



• The form of this Lagrangian was first written down by Eichten and Hill ( Eichten 



and Hill, 1990a). 



• The lattice action preserves all the continuum heavy quark symmetries discussed 
in the previous section. 
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2.2 Propagator 



From the Lagrangian eq. (2.3 1 we have the defining equation for the propagator 

: , 1 , (Vg + 6m)G h (x, y) = 6(x - y)P+ = a" 4 TT 6x, P+ . (2.5) 
1 + a dm XJ - a a 

Obviously Gh(x,y) is proportional to S(x — y). Writing Gh(x,y) = g(n ,ko;xL)5(x — 
y)P+ with Xq = ariQ , yo — afc , the above equation yields a simple recursion for 
g (no + 1, kg] x) in terms of g(no, &o; x ) which is solved by 

5(710, fc ; x) = 6(n - fc )(l + aM" ( "°" fco) ^(y, a?5 0) f , (2.6) 

V(x,x;Q) = l, V(x,y + ad;0)^V(x,y;0)U(y,0), (2.7) 

where 

tf(no-*b) = (j n ° (2-8) 

The static propagator reads 

G h (x, y) = 8(x Q - y ) <5(x - y) exp ( - 6m (x - yo)) V(y, x; 0) f P+ , (2.9) 
5m = i ln(l + aSm) . (2.10) 

The object V{x, y \ 0) parallel transports fields in the fundamental representation from 
y to x along a time-like path. Note that the derivation fixes 9(0) = 1 for the lattice 
^-function. As in the continuum, the mass counter term Sm just yields an energy shift; 
now, on the lattice, the shift is 



77,QCD j^stat 



+ m barc , m hm . c = 5m + m. (2-11) 

Sm— 



It is valid for all energies of states with a single heavy quark or anti-quark. As in the 
continuum the split between Sm and the finite m is convention dependent. 
In complete analogy the anti-quark propagator is given by 

Gy,(x, y) = % - xo) <5(x - y) exp ( - 6m (y - x )) P(x, y; 0) P_ . (2.12) 

2.3 Symmetries 

All HQET symmetries are preserved on the lattice, in particular the U(2F) spin-flavor 
symmetry and the local flavor-number conservation. The symmetry transformations 
can literally be carried over from the continuum, e.g. eq. (1.44). One just replaces the 
continuum fields by the lattice ones. 

Note that these HQET symmetries are defined in terms of transformations of the 
heavy quark fields while the light quark fields do not change (unlike e.g. standard 
parity). Integrating out just the quark fields in the path integral while leaving the 
integral over the gauge fields, they thus yield identities for the integrand or one may 
say for "correlation functions in any fixed gauge background field" . 
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2.4 Symanzik analysis of cutoff effects 

According to the — by now well tested^ — Symanzik conjecture, the cutoff effects of 



a lattice theory can be described in terms of an effective continuum theory ( jSymanzik, 
1983a Symanzik, 1983& Luscher et al., 1996). Once the terms in Symanzik's effective 



Lagrangian are known, the cutoff effects can be canceled by adding terms of the same 
form to the lattice action, resulting in an improved action. 



For a static quark, Symanzik's effective action is (Kurth and Sommer, 2001) 



S c ff = So 



Si = / d*x&(x) 



(2.13) 



where -2o(a;) — Sf^ (x) is the continuum static Lagrangian of eq. (1.47) and 



&iix) 



i=3 



CiOi(x) , 



(2.14) 



is given in terms of local fields with mass dimension [Oj(a;)] = 5. Their coefficients 
Cj are functions of the bare gauge coupling. Assuming for simplicity mass-degenerate 
light quarks with a mass mi, the set of possible dimension five fields, which share the 
symmetries of the lattice theory, is 



C>3 = ip h D D t() h , 04 = Tfi\ ip h D iph , 5 = m\ ^ h -0h 



(2.15) 



Note that P + a jP + = means there is no term ^^o-QjF^ijju, and ^j^DjDjtl)^ can't 
occur because it violates the local phase invariance eq. (1.44). Finally iph°~jkF'jk'4'h is 



not invariant under the spin rotations eq. (1.43) 



Furthermore, we are only interested in on-shell correlation functions and energies 
For this class of observables O3 



O4 do not contribute ( Luscher and Weisz, 
Luscher et al., 19961 because they vanish by the equation of motiorj^J 



1985 



D Q ip h = . 



(2.16) 



The only remaining term, O5, induces a redefinition of the mass counter-term 5m 
which therefore depends explicitly on the light quark mass. 

We note that for almost all applications, 5m is explicitly canceled in the relation 
between physical observables and one thus has automatic on-shell 0(a) improvement 
for the static action. No parameter has to be tuned to guarantee this property. Still, 
the improvement of matrix elements and correlation functions requires to also consider 
composite fields in the effective theory. 



See Peter Weisz' lectu res for a theoretical discussion and chap ter I of (Sommer, 20061 for an 
overview of tests. Finally (IBalog et al., 2009b Balog et al., 2009a) represents the most advanced 
understanding of the subject. 

2 The equations of motion follow just from a change of variable in the path in tegral. Contact terms 
are re-absorbed into the free coefficients c;. We refer to Peter Weisz' lectures or l |Luscher et al., 1996[ ) 
for a more detailed discussion. 



16 Lattice formulation 



Exercise 2.1 The static quark anti-quark potential. 
A (time-local) field 

0(t, x, y) = tp h (x) V(x, y)-y 5 ifc(y) , x = yo = t 

with V(x, y) being a parallel transporter from y to x in xo = t plane, can be used to annihilate 
a quark-anti-quark pair at a separation x — y, while 



0(t, x, y) = -V>h(f) P(j/> x )7s V>h(» , x = y = i 

will create a quark-anti-quark pair at a separation x — y. 
Show that for t > 



(2.17) 



( 0(t, x, y) O(0, x, y) ) = const. e" 2t 4m VK(i, x - y) 



(2.18) 



where W is the Wilson loop introduced in the lectures of P. Hernandez. Since the energy 
levels of HQET are finite (after inclusion of a suitable 5m), one can conclude that 



Vk(x-y) 



lim d t ln(W(t, x - y)) + 25m 



(2.19) 



is a finite quantity: the divergent constant in the bare potential is absorbed by 8m, i.e. by a 

renormalization of the heavy quark mass. 

Furthermore, from the 0(a) improvement of HQET, one concludes (Nccc cTand Sommer,| 



2002) 



V R (x - y) = Vg»*(r = |x - y|) + 0(a 2 ) 
if the action for the light fields is O(a) improved. 



(2.20) 



2.4.1 Renormalized and improved axial current. 

We now also have to specify the discretization of the light quark field ip. We will 



generically think of a standard 0(a)-improved Wilson discretization ( Sheikholeslami 
and Wohlert, 1985]|Luscher et al, 1996 ) but occasionally mention changes which occur 



1998 


Liischer, 1998 


3 


or a Wilson rej 


jularization with a twisted 


et al., 2001a 


Frezzotti et al, 20016 


Frezzotti and Rossi, 2004 



study the time component of the axial current. In Symanzik's effective theory it is 
represented by 



k=l 



stat 
)k 



) k , < at = ^ 7 o75^h 



(2.21) 



with some coefficients cu^. Here the flavor index of the field ip is suppressed. It is 
considered to have some fixed but arbitrary value for our discussion, except where we 
indicate this explicitly. A basis for the dimension four fields {{SAQ tat )k} is 



3 When the chiral symmetry realization of domain wall fermions (Shamir, 1993) is good enough, 
these fermions can of course also be considered to have an in practice exact chiral symmetry. 
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(SA s tat ) 1 = V^7,75^h , (SAf a % = V>75A^h , 

(2.22) 

(<L4 s tat ) 3 = ^ o75 ^h , (SA s ta % = my ^ 7 o75^h . 

From eq. ( |2. 16 ) we see that k — 2 does not contribute, while the equation of motion 
for if; relates ((L4o tat ) 3 , (<L4o tat ) 4 and (SAf at )i. We choose to remain with k = 1 (and 
in principle k = 4), but for simplicity assum^jami -C 1; we can then drop (8A * at )4. 
So for on-shell quantities the effective theory representation is 

« at )cff = Af at + awit&C*)! . (2.23) 

In order to achieve a cancellation of the 0(a) lattice spacing effects, we add a cor- 
responding combination of correction terms to the axial current in the lattice theory 
and write the improved and renormalized current in the form 

(i4g»*)o = Zf at (g , an) (Af at ) , (2.24) 
(Af at ) = A#* + acf* (g ) ^ 75 §(^ + V*)V„ , (2.25) 

with a mass-independent renormalization constant Z^ at and a dimensionless improve- 
ment coefficient, cj at , depending again on go but not on the light quark mass. 

The improvement coefficients can be determined such that for this (time component 
of the) improved axial current we have the representation 

« at )cff = ^TtoTStfh + 0(a 2 ) , ( 2 - 2 6) 

in the Symanzik effective theory. In other words u)\ is then 0(a) and cutoff effects are 
0(a 2 ). 

The symmetries of the static theory are strong enough to improve all components 
of the flavor currents in terms of just c^ at and to renormalize them by Z^ at . Let us 
discuss how this works. 



2.5 The full set of flavor currents 

The previous discussion literally carries over to the time component of the vector 
current, 

V stat = ^7oV>h ■ (2.27) 

Its improved and renormalized lattice version may be chosen as 

(Vjf at ) = Z^ tat (Vf tllt )o (2.28) 
(Vf tat ) = Vn^ h + acT^li \ (Vj + V*) ■ (2.29) 

4 If the light quark action has an exact chiral symmetry or the light quarks are discretized with 
a twisted mass term at full twist, this restriction is unnecessary, since the term is excluded by the 
symmetry. Note that (<5Ag tat )i is, however, not forbidden by chiral symmetry and c^ at is necessary 
for 0(a)-improvement in any case. 
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The chiral symmetry of the continuum limit can be used to relate Zy at , Cy at to 
^Btat ^ c stat j n ^ ne f ii owm g wav . We assume A/f > 2 massless light quarks. Then the 
infinitesimal transformation 

SH>(x) = ir^sVW, W^) = Vfahs^ , (2.30) 

with the Pauli matrices r a acting on two of the flavor components of the light quark 
fields ip, ip, is a (non-anomalous) symmetry of the theory. Identifying V^ stat = V'i7oV'h> 
where ip 1 is the first flavor component of ip, the vector current transforms as V^ stat = 
— |AQ tat . The same property can then be required for the renormalized and improved 
lattice fields, 

6%(Vg«) = -|« at )o + 0(a 2 ) . (2.31) 

This condition can be implemented in the form of Ward identities relating different 
correlation functions, in particular in the Schrodinger functional . We refer to A. 



Vladikas' lectures and (Liischer, 1998) for the principle; practical implementations 



have been studied in (Hashimoto et al, 2002 Palombi, 2008). Such Ward identities 
determine Zy at , Cy at in terms of Zf at , d£ at 



Furthermore by a finite spin-symmetry transformation (with uu of eq. (L4TJ) 

ih -> = e- i7rCT */ 2 Vh = -io- k ^ h , rf b = , (2-32) 

we have 

and we can require the same for the correction terms, 

[«5y stat ] ' = SAf at , [SAf at ] ' = SV£ tat . (2.34) 

We leave it as an exercise to determine the form of 8A s £ at , 8V^ tat . The discussed 
transformations are valid for the bare lattice fields at any lattice spacing. Thus renor- 
malization and improvement of the spatial components is given completely in terms of 
the time-components once we define the renormalized fields to transform in the same 
way as the bare fields. A last property to note before writing down the renormalized 
and improved fields is that we have 

Z^(g , an) = Zy=f A (g ) Zf*\g , afx) (2.35) 

with a /^-independent function ZyJ^go) and up to 0(a 2 ), as soon as we require 
eq. (pljl.pl 



5 A formal argument is as follows. Rewrite eq. 1 2.31 1 in terms of the bare operators, 
Zy^igo, a i l )^iyS taX = -|- Z A at (S'0,aA t )^-o tat + 0(a 2 ). Since the bare, regularized, operators 
yjstat ^stat carr y no ^-dependence, we see that Zy at (30, afi)/Z^ at (go, iaO is a function of 30 only, 
apart from 0(a 2 ) cutoff effects. To make the argument more rigorous one should rewrit e the equation 
in the form of correlation functions which represent a Ward identity equivalent to eq. ( |2.31| . 
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Let us disregard the 0(a) improvement terms for simplicity. We can then summa- 
rize what we have learnt about the renormalization of the static-light bilinears as 

(A s R tat )o = Z^(g , a?) A s tat , (2.36) 

(V£ tat )o = ZI tat ( 3o , ail ) Z$f A (g ) y stat , (2.37) 

(VS tat ) fe = Zl tat ( ff0 ,a/i)yf at , (2.38) 

{A^% = ZT\ 90 , a M ) Z$f A (g ) Af at , (2.39) 



where Zyjj^(gp) can be determined from a chiral Ward identity (Hashimoto et al, 



2002||Palombi, 2008[ ). Note that we denote the flavor currents in HQET in complete 
analogy to QCD. Still they do not form 4- vectors, as 4-dimensional rotation invariance 
is broken in HQET. For example (Ap tat ) cannot be rotated into (A^ 3jt ) k by a 90 degree 
lattice rotation. 

The only bilinears which are missing here are scalar, pseudo-scalar densities (and 
the tensor). These are equivalent to A^ at and V' stat in static approximation, for ex- 
ample 

HsA = ^757o^h = -< at , = ^rotfh = Kf at • (2.40) 

At this stage it is therefore unnecessary to introduce renormalized scalar and pseudo- 
scalar densities. 

We have so far written down expressions for the relevant renormalized heavy-light 
quark bilinears. The Z-factors can be chosen such that correlation functions of these 
fields have a continuum limit (with 5m, gauge coupling and light quark masses properly 
determined). Beyond this requirement, however, also the finite parts need to be fixed 
by renormalization conditions. We have fixed some of them such that the renormalized 
fields satisfy chiral symmetry and heavy quark spin symmetry. Only one finite part 
(in Z s £ at ) then remains free. Preserving these symmetries by the renormalization is 



natural, but not absolutely required; e.g. eq. (2.33) could be violated in terms of the 
renormalized fields. As long as one just remains inside the effective field theory these 
ambiguities are not fixed. The proper conditions for the finite parts, valid for HQET 
as an effective theory of QCD, have to be determined from QCD with finite heavy 
quark masses. We will return to this later. 

We may, however, already note that for renormalization group invariant fields, 
these ambiguities are not present. The renormalization group invariants are thus very 
appropriate. Still, relating the bare lattice fields to the renormalization group invari- 



ant ones is a non-trivial task in practice (Liischer et al, 1991 Capitani et al, 1999). 



We will briefly discuss how it can be done (and has been done) for the static-light 
bilinears ( |Kurth and Sommer, 20"0l]|Heitger et al, 2003||Della Morte et al, 20076) ). For 



this and other purposes we need the Schrodinger functional . In the following we just 
give a simplified review of it and describe how static quarks are incorporated. Some 
more details are discussed by Peter Weisz. 

2.6 HQET and Schrodinger Functional 



The Schrodinger functional ( |Symanzik, 1981||Luscher et al, 1992|Sint, 1994jSint, 1995} 
can just be seen as QCD in a finite Euclidean space-time of size T x L , with spe- 
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cific boundary conditions. It is useful as a renormalizable probe of QCD, providing a 
definition of correlation functions which are accessible at all distances, short or long: 
gauge invariance is manifest and even at short distances (large momenta) cutoff effects 
can be kept small. It will help us to perform the non-perturbative renormalization of 
HQET and its matching to QCD. In all these applications it is advantageous to have 
a variety of kinematics at ones disposal. One element is to have access to finite but 
small momenta of the quarks (think of the free theory, a relevant starting point for 
the short distance regime). 

To this end, the spatial bounda ry conditions were chos en to be i])(x + Lk) = 



ip(x) , ip(x + Lk) = e k ip(x) m (Sint and Sommer, 1996), which allows momenta 



pk = *nh + 0k_ > he7L ^ (241) 

in particular small ones when Ik = 0. Performing a variable transformation ip(x) —> 
e l6hXk / L ip(x), ip(x) — > e~ l9kXk / L ip(x), for < Xk < L — a, we see that this boundary 
condition is equivalent to periodic boundary conditions (without a phase) for the new 
fields, while the spatial covariant derivatives contain an additional phase, for example 

V k ip(x) = - [e Wka / L U{x, n)ip(x + ak) - tjj(x)] , (2.42) 



see also Sect. A.l The phase 8ka/L can be seen as a constant abelian gauge potential 
and the above variable transformation as a gauge transformation. Of course, the angles 
8k which we will set all equal from now on (8k — 8) , are not specific to the Schrodinger 
functional ; they just have first been used in this context. 



The standard Schrodinger functional boundary conditions in time are (Sint, 1994 



Sint, 1995) 



P + V(aOL=o=0, P_V(aOU 0=T = 0, (2-43) 

and 

^(x)P. | xo=0 = , i>{x)P + \ X0 = T = . (2.44) 

The gauge fields are taken periodic in space and the space components of the contin- 
uum gauge fields are set to zero at xq — and xq — T (on the lattice the boundary 
links U(x,k) are set to unity)j^] 

For the static quark the components projected by P_ vanish anyway, so there is 
just 

P+^(aO| Xo =o=Q, ^ h (z)P+|x o =T = 0. (2.45) 

Defining 

ip h (x) = if xo < or x > T, (2.46) 

the lattice action for the static quark with Schrodinger functional boundary conditions 
can be written as 



In I Luscher et al., 1992[|Liischer et al., 1994 1 the definition of a renormalized coupling uses more 



general boundary conditions tor the gauge fields, but these are not needed here 
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S h = 1 - X — a 4 V'h (») [ Vo + M (a) (2.47) 



as before. In general the improvement of the Schrodinger functional requires to add 
boundary terms to the action as a straightforward generalization of Symanzik improve- 
ment. These terms are dimension four composite fields located on or at the boundaries, 



summed over space (Liischer et ai, 1992 Liischer et al, 1996). Since they are not so 



important here and are also known sufficiently well, we do not discuss them. We just 



note that no boundary improvement terms involving static fields are needed (Kurth 



and Sommer, 2001), since the dimension four fields vanish either due to the equation 
of motion or the heavy quark symmetries. 

We take the same periodicity in space as for relativistic quarks, 



ip^x + Lk) = iph(x) , ip h (x + Lk) = ip h (x) . 



In the static theory this has no effect, since quarks at different x arc not coupled, but 
it plays a role at order 1/m where is a useful kinematical variable. 

An important feature of the Schrodinger functional is that one can form gauge 
invariant correlation functions of boundary quark fields. In particular, one can project 
those quark fields to small spatial momentum, e.g. p = 1/L X (6, 6, 0) for the quarks 
and — p for the anti-quarks. For the precise definition of the boundary quark fields we 
refer to (Liischer et ai, 1996 ) or for an alternative view we refer to ( |Liischer~0 06 ) . The 



details are here not so important. We only need to know that these boundary fields, 
i.e. fermion fields localized at the boundaries, exist. Those at xq = are denoted by 



Ci(x), d(x), Ch(x), Ch( X ): 

and those at xq = T by 

Ci'(x), Ci'W, Ch'(x), Ch'(x). 
2.6.1 Renormalization 

These boundary quark fields are multiplicatively renormalized with factors Zq , Z^ h , 
such that (Ci(x))r = Z^Ci(x) etc. 

To illustrate a first use of the Schrodinger functional and the boundary fields we 
introduce three correlation functions 



22 Lattice formulation 



/l tat (*o,0) 



a 6 x - 



((Af at ) ( 2; )C h (y)75Ci(z)> 



(2.49) 



/f at W = -|p E <C 1 , (u)7 5 Ch , (v)C h (y)75Ci(^)) 



(2.50) 



a 

"2L 5 



xi,x 2 ,y,z 



(Ch'(x) 7 5Ch'(o) C h (y)7 5 Ch(z)) 



(2.51) 



In the graphs, double lines are static quark propagators. Note that the sum in eq. ( 2.51 ) 



runs on X\ and X2 and therefore yields an cc 3 -dependent correlation function. We 
further point out that ^ y etc. project the boundary quark fields onto zero (space) 
momentum, but together with the abelian gauge field, this is equivalent to a physical 
momentum pk = 9/L. For example the time-decay of a free mass- less quark prop agator 
projected this way contains an energy E(9/L,6/L,6/L) = y/36/L, cf. eq. Q >. 
The above functions are renormalized as 



r .fStatl VStat ry ry /*stat [PStatl ryl ryl rstat X -fllhl ry\ J"hh 

L/A Jr = z a ^Ch^C /a » [h J R - z Ch Z CJi ' LAJ R -^Ch/i 



(2.52) 



We remind the reader that an additional renormalization is the mass counter-term of 
the static action. 
The ratio 



/f*(r/2,0) 



(2.53) 



renormalizes in a simple way and also needs no knowledge of 5m, since it cancels out 



due to eq. (2.9). It is hence an attractive possibility to define the renormalization 



constant through this ratio. Explicitly we may choose 



/l tat (W) 



fi=l/L, T = L. 



1=0 



(2.54) 



which defines the finite part of Z^ at in a so-called Schrbdinger functional scheme. As 



usual the factor 



f^(L/2,6)/^f? 



3=0 



is inserted to ensure Z^ tat 



l + O( ff0 2 ). 



The name Schrodinger functional scheme just refers to the fact that the renormal- 
ization factor is defined in terms of correlation functions with Schrodinger functional 
boundary conditions. While Z^ at refers to a specific regularization, the renormal- 
ization scheme is independent of that and can in principle be applied in a continuum 
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regularization. Many similar Schrodinger functional schemes can be defined (e.g. (Heft 
ger et al, 2003)), but by the choice T = L, 9 = 0.5 we have made eq. (2.54) unique. 
It is implied that the light quark masses are set to zero. We will soon come back to 
the /^-dependence of the renormalized current and its relation to the RGI current. 
First let us show some numerical results which provide a non-perturbative test of the 
renormalizability of the static theory. 



2.7 Numerical test of the renormalizability 

The above listed renormalization structure of the Schrodinger functional correlation 
functions is just deduced from a simple dimensional analysis. A number of 1-loop 



calculations of the correlation functions defined above as well as of others ( 


Kurth 


and Sommer, 2001|Kurth and Sommer, 2002 


Delia Morte et al, 2005 


Palombi, 2008 1 


confirm the structure eq. (2.52 


) and more generally the renormalizability of the theory 



(by local counter-terms). 



Also non-perturbative tests exist. A stringent and precise one (Delia Morte et al, 



2005 ) is based on the ratios 



Ha(0,0') 



/i tat (T/2,0) 



/I tat (r/2, 



6 



/f at W 



h{d/L,6) 



fr(d,9) 



(2.55) 

The additional dependence on L and the lattice resolution a/L of these ratios is not 
indicated explicitly. With eq. (2.52), we see that all renormalization factors cancel 

0, approached asymptotically 
where L is kept fixed in units of the 



in these ratios. They should have a finite limit a/L 
with a rate (a/L) 2 . This is tested in Fig. 



2.1 



reference length scale ro(Sommer, 1994) to L/r$ = 1.436. This choice corresponds to 
about L 0.7 fm. The same continuum limit has to be reached for different lattice dis- 
cretizations. Also this universality is tested in the graphs, where four different choices 
of the covariant derivative Dq in the static action are used. All actions defined by the 
different choices of Dq have the symmetries discussed earlier. 



2.8 Scale dependence of the axial current and the RGI current 

Let us first recapitulate the scale dependence in perturbation theory. At one-loop order 
one has 



Z% at (g , (m) = l+g 2 [B - 70 ln(o (i)] + 



7o 



1 

"4^2 



(2.56) 
(2.57) 



In the lattice minimal subtraction scheme the Z-factors are polynomials in ln(a/i) 
without constant part; thus Bq = 0. Instead, when the renormalization scheme is 



defined_by eq. (|2.54|) a one- loop computation of /A tat , /f tat yields(|Kurth and Sommer, 



|2()()1 ■ p] />';; = —0.08458. As usual there is the renormalization group equation (RGE) 
(remember g 2 — g 2 + 0(g 4 )) 



7 As usual in perturbation theory, terms of order (afi) n , n > 1 are dropped. 
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0.005 



0.015 (a/L) 2 



Fig. 2.1 Lattice spacing dependence of various ratios of correlation functions for which Z-factors 
cancel. Different symbols correspond to different actions. Computation and figure from | |Della Morte| 
\et gl, 2005] |. 

M^« at )o = 7(5)(^ at )o , l(g) = -f {lo+fli + ...}■ (2.58) 



Combining it with the RGE for the coupling eq. (A.29) it is easily integrated to (see 
eq. (A. 32) for the definition of the beta-function coefficients hi) 



(OoM = (^ RGI )o exp 



dec- 



(A RGI ) [26 3 2 ] 7o/2bo 



0(x) 



exp 



da; 



(3(x) box 



(2.59) 
(2.60) 



where eq. (2.60) provides the definition of the lax notation for the second factor in 



can also be written as 



eq. (2.59). The integration "constant" is the renormalization group invariant field. It 



(A RGI ) - lim [2b g 2 (/,)}' 

fl — ^oo 



since the last factor in eq. ( 2.60 1 converges to one as /i — > oo. Using also that 70, 60 are 
independent of the renormalization scheme, as well as Os(/i) = Os'(//)(l + 0(ff 2 (A*)) 
(valid for any operator O and standard schemes S, 5"), this representation also shows 
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that the renormalization group invariant operator (A )o is independent of scale and 
scheme 

Let us now go beyond perturbation theory and start from a non-perturbative defi- 
nition of the renormalized current, such as eq. (2.54), together with a non-perturbative 
definition of a renormalized coupling ( Liischer et al, 1991|Liischer et al, 1992|Liischer 



et al, 1994). With the step scaling method discussed in more detail by Peter Weisz, 



one can then determine the change 

(Ar t )o(^) = < tat (5 2 (2M))(^ tat )o(2/i), li = l/L (2.61) 

of the renormalized field (A^ at )o(n) when the renormalization scale [i is changed by a 
factor of two. The so-called step scaling function a^ at is parameterized in terms of the 
running coupling g(fJ-). Its argument is fi = l/L in terms of the linear extent, L = T, 
of a Schrodinger functional . 

Instead of the scale dependence of (A^ lit )o(n) we will often discuss a generic matrix 
element 

$00 = (a\(^) (ji)m (2.62) 

of the associated operator A^ at in Hilbert space. 

In a non-perturbative calculation, the continuum cr^ at is obtained through a nu- 
merical extrapolation 

af &t (u) = Urn E s A tat (u,a/L) (2.63) 

a/L— >0 



of the lattice step scaling functions 

Z*ta.t(„ n „lo.T\ 

(2.64) 



Er(, a /D- zr(s »' 8/2£) 



Zf at (g ,a/L) 



g 2 (i/L)=u 



obtained directly from simulations. Here, g 2 (l/L) is kept fixed to remain at constant 
L while L/a is varied in the continuum extrapolation. 

The /i-dependence of $ can then be constructed iteratively via 

u Q = g 2 (l/L ) , L n =2 n L , : $(l/L„+i) = <4 tat (u n ) $(1/L„) 

u n +i = cr(u n ) , 

where the step scaling function a of the running coupling enters. The length scale Lq 
is chosen deep in the perturbative domain, typically Lq w 1/100 GeV and therefore 
the ^(-dependence can be completed perturbatively to infinite fj,, i.e. to the RGI using 
eq. ( |2~60l ). 

For Nf = 2 the analysis has been done for /j, w 300 Me V... 80 GeV. After it was 
verified that the steps at smallest L (L < L%) are accurately described by perturbation 
theory (see Fig. 2.2), the two-loop anomalous dimension was used in eq. (2.60) with 



Of course, a trivial definition dependence due to the choice of pre-factors in eq. I 2.60 i is present. 
Unfortunately there is no uniform choice for those in the literature. 
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Fig. 2.2 Relation 4 , (At)/4 , RGl between RGI and matrix element at finite /i in a Schrodinger func- 
tional scheme and for Nf = 2. The A-parameter in the SF-schemc is around 100 MeV. Everything 
was computed non-perturbatively from continuum extrapolated step scaling functions ( |Della Morte| 
\et al, 2007b} . 



[i = \jL v to connect to the RGI current. The result ( Kurth and Sommer, 2001[[Heitger 
et al., 2003 Delia Morte et al., 2007 b) is conveniently written as 



^A,RGl(5o) 



RGI 



xZt^(g , ai i) 



(2.65) 



p=l/(2L max ) 



where only the second factor depends on the lattice action, and g 2 (l/£ max ) = it max is 
a convenient value covered by the non-perturbative results for the above recursion. 

We show the result for the first factor in Fig. |2.2| for a series of /i with Nf = 2 
dynamical quarks. The different points in the graph correspond to different n in the 
recursion. Note that the two- loop running becomes accurate only at rather small L. 
There is an about 5% difference in between a two-loop result and the non- 

perturbative one at the smallest fi. 

The details of this calculation and strategy is not that important for the following. 
We have mainly discussed it since 

- first the RGI matrix elements play a prominent role in HQET in static approxi- 
mation and it is relevant to understand that they can be obtained completely non- 
perturbatively and 

- second we also want to later emphasize the difference between the here used - by 
now more classic - renormalization of the static theory and the strategy discussed in 
Sect. El 

Let us further note that the strategy above has been extended to four-fermion 
operators relevant for B — B oscillations in (Palombi et al., 2006 Dimopoulos et al., 
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Fig. 2.3 Relation c = ^rgi/^Cm) of the RGI matrix clement 'I'rgi and the matrix element at 
finite fi of the "V A + AV" four— fermion operator in a Schrodinger functional scheme. It was computed 
non-perturbatively from continuum extrapolated step scaling functions l |DimopouIos et ql., 2008} . 



2008 ). Since more than one operator is involved and twisted mass QCD is used in order 
to avoid mixing with operators of wrong chirality, the strategy and computation are 
somewhat more involved. In particular two different four-fermion operators contribute 
at leading order in I/m. As an example, we just show in Fig. 2.3 the result for the 
four-fermion operator which dominates in the physical process. 

2.9 Eigen-states of the Hamiltonian 

The eigenstates of the static Hamiltonian can be diagonalized simultaneously with the 
local heavy flavor number operator (remember Qh(x) — ip h (x)iph(x) )• We consider 
a finite volume with periodic boundary conditions. Since the theory is translation 
invariant, there is a k x L 3 /a 3 - fold degeneracy of states with a single heavy quark|^J 
where k arises from degeneracies on top of the translation invariance discussed here. 
For the lowest energy level one can choose a basis of eigenstates of the Hamiltonian as 

|B(x)), (B(x)|S(y))=2<S(x-y), Q h (y)|S(x)) = *(x - y) (2.66) 

or their Fourier transformed 

|fl(p)) = a 3 ]T e-^l^x)) , (B(p')\B(p)) = 2(2n) 3 6(p - p') (2.67) 



5(p-p') = (L/(2n)) 3 l[5 lil ,, h 



2-kU 

17 



G TL. 



(Here we set 9 = 0.) We usually work with the zero momentum eigenstate, denoted 
for short by \B) = \B(p — 0)), as this is related to an eigenstate of the finite mass 
QCD Hamiltonian, which in finite volume has normalization (B\B) = 2L 3 . 



3 In certain types of quark smearing, this has to be properly taken into accounti Christ et al, 2007) 



3 

Mass dependence at leading order in 
1/m: Matching 



We now discuss the "matching" of HQET to QCD using the example of a simple 
correlation function. As mentioned before, the issue is to fix the finite parts of renor- 
malization constants such that the effective theory describes the underlying theory 
QCD. Throughout this section we remain in static approximation. Matching including 
1/m terms will be discussed in the next section. 

3.1 A correlation function in QCD 

We start from a simple QCD correlation function, which we write down in the lattice 
regularization, 

Cg£(*o) (4>(*)4(0)) (3.1) 

X 

with the bare heavy-light axial current in QCD, — ipj^lb^b, and A* = V'bT^TsV'- 
The current is formed with the relativistic b-quark field ip^. In QCD, the renormaliza- 



tion factor, Zx{go), is fixed by chiral Ward identities (Bochicchio et al, 1985 Luscher 



et ai, 1997). It therefore does not depend on a renormalization scale. 



One reson to consider this correlation function is that at large time the B-meson 
state dominates its spectral representation via 

CgW = ZWY,ml(x)\B)^(B\A a (0)\0)e-^ [l + (e^ A )] 

X 

= Z|i(0|4(0)|B)(i?|Ao(0)|0)e-^" lB [l + 0(e-« A )] (3.2) 

and the B-meson mass and its decay constant can be obtained frornP] 

TaS D (^o) = -do ln(C^ D (z )) = m B + 0(e~*° A ) (3.3) 
[d>Q CD ] 2 EE/Im B (3.4) 



(B\Z A A o \0) =2 lim exp(i r^(;Eo))CS£ u (a:o) 

Xq— >00 



Note that we use the normalization eq. (1.52) (in finite volume (B\B) = 2L 3 ) for the 
zero momentum state \B). The gap A is the energy difference between the second 

1 It is technically of advantage to consider so-called smeared-smeared and local-smeared correlation 
functions, but this is irrelevant in the present discussion. 
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energy level and the first energy level in the zero momentum (flavored) sector of the 
Hilbert space of the finite volume lattice theory. 

3.2 The correlation function in static approximation 

In the static approximation we replace Z&A — > Z^ a,t (g , fia)AQ tat and define 

Cft R (x ) = (zr^Cf^ixo) = (zr') 2 « 3 E (< at (z)« at ) f (0)) (3.5) 

— \ / stat 

X 

T s lf(xo) = -doHCfl t (x )), (3.6) 
[$( M )] 2 = |(B|Zl tat < at |0) stat | 2 (3.7) 
= 2 lim exp^oT^^o))^) 2 ^^). 

Xq— fOO 

The /x-dependence of $ results from the renormalization of the current in the effective 
theory, 

ZX tat ( 5o , M a) = 1 + g 2 a [B - 7o ln(a M )] + O( 5o 4 ) . (3.8) 

Different renormalization schemes have different constants Bq. Alternatively one uses 
the renormalization group invariant operator (A I )o- We come to that shortly. 

3.3 Matching 

The correlation function C^ A D and the matrix element $Q CD , eq. (3.4), are indepen- 



dent of any renormalization scale, due to the chiral symmetry of QCD in the massless 
limit. But of course they depend on the mass of the b-quark. 

In the effective theory we first renormalize in an arbitrary scheme, which we do not 
need to specify for the following, resulting in a scale-dependent <I>(/i). The two quan- 
tities arc then related through the matching equation (without explicit superscripts 
"QCD" we refer to HQET quantities, here static), 

$Q CD (m) = C match (™, /i) x + 0(l/m) . (3.9) 

Somewhat symbolically the same equation could be written for the current instead of 
its matrix element; we write "symbolically" since the two currents belong to theories 
with different field contents. However, thinking in terms of the currents, it is clear 



that eq. (3.9) can be thought of as a change of renormalization scheme in the effective 



theory, where the new renormalization scale is m = and the finite part is exactly 



fixed by eq. (3.9). In fact, since at tree level we have constructed the effective theory 
such that $ = <I>Q CD , the tree-level value for C matc h is one and we have a perturbative 
expansion 

C mat ch("i, A*) = 1 + c i( TO /M)ff 2 (/ i ) + • • ■ (3.10) 



The finite renormalization factor C matc h may be determined such that eq. (3.9) holds 
for some particular matrix element of the current and will then be valid for all matrix 
elements or correlation functions. Some aspects of the above equation still need ex- 
planation. The /i-dependence is not present on the left-hand-side and this should be 
made explicit also on the right-hand-side; further one may wonder which definition of 
the quark mass and coupling constant one is to choose. 
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3.3.1 One-loop 



Before coming to these issues, it is illustrative to write down explicitly what eq. (3.9) 



looks like at 1-loop order. Ignore for now how we renormalized the current in eq. (2.54) 
and use instead lattice minimal subtraction, 

ZI tat ( 5o , M a) = 1 - 7 oln(a^).9o + ••• , (3.11) 

in the static theory. 

Instead of the decay constant, use as an observable a perturbatively accessible 
quantity. We tak^] 

$ qcd = Y *c»^ mR L) s B ZA{go) MW' m f > (3-12) 

where for our one-loop discussion we do not need to specify the normalization condition 
for the renormalized heavy quark mass fn = tor and coupling g = gn- The one- loop 
expansion of these functions has been computed ( Kurth and Sommer, 2002] ), and the 
result can be summarized as 



$ = Y^(6,fx,L) = hmZt at (g ,H ^S=ff (3.13) 
= A{6)[\ - 70 HliL) gl] + D(0)<& + 0(4) 

in static approximation and 

0>QCD = A(6)[l + (£>' - 7o ln(TO R L))^J + D{6 )gl + 0(l/(m K i)) + O(^) 

in QCD. From these expressions we can read off 

ci (tjir/m) = 7o In(^/m R ) + D' . (3.14) 

Furthermore, the fact that the same functions A(8), D{&) appear in the static theory 
and in QCD is a (partial) confirmation that the static approximation is the effective 
theory for QCD. In particular the logarithmic L-dependence in QCD matches the one 



in the static theory. With eq. (3.14), the matching of QCD and static theory holds for 



all 9, and also for other matrix elements of A . 
3.3.2 Renormalization group invariants 

Having seen how QCD and effective theory match at one-loop order, we now proceed 



to a general discussion of eq. (3.9), beyond one-loop. Obviously, the /^-dependence in 



eq. (3.9) is artificial, since we have a scale-independent quantity in QCD. Only the 



2 The correlation functions /a, fi are the relativistic versions of / A tat , /| tat . 
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mass-dependence is for real. We may then choose any value for /i. For convenience we 
set all renormalization scales equal to the mass itself 

H = m* = m(m*) , = ff(m*) , (3.15) 

where ro(^i), <?(u) are running mass and coupling in an unspecified massless renormal- 
ization schemer] This simplifies the matching function to 

C ma tch(m^, m*) = C ma tch(ff*) = 1 + ci(l) gl + ... . (3.16) 

Further we want to eliminate the dependence on the renormalization scheme for 
fn,g, (A^ at )o. As a first step we change from <&(//) to the RGI matrix element 

( raW <y( x ) 1 
$ RG i = exp <^ - / dx—} $(//) , (3.17) 



and arrive at the form 



$ QCD = C ma tch(<7*) x $( M ) = CUtch^*) exp J /"* dx^- } $ RG i (3.18) 

[J /3(a;) J 

f 7match(a;) 1 

sexp {/ di ^w-r Eo " <319) 

Everywhere terms of order 1/m are dropped, since we are working to static order. 
Eq. (3.19) defines 7 ma tch, which describes the physical mass dependence via, 

^QCD ^ = 7match(ff*) , (3.20) 

but it still depends on the chosen renormalization scheme through the choice of m (the 
scheme, not the scale). We eliminate also this scheme dependence by switching to the 
RGI mass, M, and the A-parameter, 

A f f s{tl) 1 ] 

exp {- / da; } , (3.21) 



expl- dx — } . (3.22) 



m (M) [ J P(x) 

Exact expressions, defining the constant parts in these equation, are given in the 
appendix. 

3 Note that is implicitly defined through m* = m(m*). 

4 In a massless renormalization scheme, the renormalization factors do not depend on the masses. 
Consequently the renormalization group functions do not depend on the masses. 
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Just based on dimensional analysis, we expect a relation 

$Q CD -Cp S (M/A)x$ RG i (3.23) 

to hold. Indeed, remembering eq. (3.15), fi = = fn, we can combine eq. ( 3.21[ ) and 
eq. ( |3.22[ ) to 

A f rs*(M/h) 1-t(x)) 

T-=cxpt- dx — } , (3.24) 

from which g+ can be determined for any value of M/A; we write = g+(M/A). It 
follows that 

dM 1 - r(<7*) V 1 

and the matching function is 

C PS (M/A) = exp I J dx v | . (3.26) 
We note that the dependence on M is described by a function^] 

= 7^atch(S*) (M/A) . (3 . 31) 

A 1 - T (5*) 



M <93> 



m acp S 

A ~ Cps^A^ 



With 



7match(3*) ° - Ivgl " 7i iatch 3* + • • • , m M ~" - bog" + ■■■ (3.32) 



we can now give the leading large mass behavior 

C PS M -°° (26 . g 2)-To/2bo „ [l n (M/A)]T°/ 26 ° . (3.33) 

Functions such as Cps convert from the static RGI matrix elements to the QCD matrix 
element; we call them conversion functions. 



5 This is seen from 

M d$_ = A^dnu m*_9$_ = 7match(g*) 
* dM m* dM $ dm* 1 - r(g*) 

i -4- , 7match(9*) 



(3.27) 



where we used 



m* = Af exp 



rs* t(x) ] 

/ da: } (3.28) 



9m* m* T (9*) dq* m± T (.9*) dm* 

= — - H — m* = — - H flfa*) *, (3.29) 

dM M p^) dM M pr g ^) vy ' dM K ' 



which shows that 



M dm* 1 



m* dM 1 - r(g+) 



(3.30) 
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Fig. 3.1 Cps estimated in perturbation theory. For B-physics we have Ajjg/Afb ~ 0.04. Figure 
from (Heitger et ai, 2004). 

An interesting application is the asymptotics of the decay constant of a heavy-light 
pseudo-scalar (e.g. B)j^] 



Fk M ^ Kapw* x[1 + 0(|ln(M/A)rl] 



(3.34) 



At leading order in 1/m the conversion function Cps contains the full (logarithmic) 
mass-dependence. The non-perturbative effective theory matrix elements, $rgIj & re 
mass independent numbers. Conversion functions such as Cps are universal for all (low 
energy) matrix elements of their associated operator. For example 



/~<QCD / \ 

^aa,r\ x o) 



x »l/m 



[C PS (4r) ^A ta RGi] 2 (< at (^) t < at (0)) +O(^) , (3-35) 
CJlH^o) (bare) 



is a straight forward generalization of eq. (3.9) 



Analogous expressions for the conversion functions are valid for the time compo- 
nent of the axial current replaced by other composite fields, for example the space 



components of the vector current. Based on the work of 


Broadhurst and Grozin, 


1991 


Shifman and Voloshin, 1987 


Politzer and Wise, 1988 


and recent efforts their 



perturbative expansion is known including the 3-loop anomalous dimension 7 ma tch ob- 
tained from the 3-loop anomalous dimension 7 (Chctyrki n and Grozin, 2003[) in the 
MS-scheme and the 2-loop matching function C mat ch ( Ji and Musolf, 1991||Broadhurst 



and Grozin, 1995 Gimenez, 1992). 



Figure 



3.1 



seems to indicate that the remaining 0(<? 6 (TOb)) errors in Cps are rel- 
atively small. However, as discussed in more detail in App. |A.2| such a conclusion is 
premature. By now ratios of conversion functions for different currents arc known to 



even one more order in perturbation theory (Bekavac et ai, 2010). We show an exam- 



ple in the first column of Fig. 3.2 where the x-axis is approximately proportional to 
gl(M/A) and for B-physics one needs l/ln(Ajjg/Mb) ~ 0.3. For a quark mass around 



6 Note the slow , logarit hmic, decrease of the corrections in eq. 1 3.34 . We will see below, in the 
discussion of Figs. [3~1|3.2| that the perturbative evaluation of Cps(M^/ A) is somewhat problematic. 
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Fig. 3.2 The ratio Cps/Cy, evaluated in the first column as described here. In columns two 
and three the expansion in is generalized to an expansion in /s), see App. |A~.2,2| The last 

column contains the conventionally used Cratch ( m Q> m Qi m Q ) /*-^match ( m Q ' m Q> m Q)i see App. A. 2 
For B-physics we have Ay^/M-^ 0.04 and 1/ \n{Ay^/ M\,) Rj 0.3. The loop order changes from 
one-loop (long-dashes) up to 4-loop (full line) anomalous dimension. 

the mass of the b-quark and lower, the higher order contributions in perturbation 
theory do not decrease significantly and perturbation theory is not trustworthy. It 
seems impossible to estimate a realistic error of the perturbative expansion. Only for 
somewhat higher masses the expansion looks reasonable. 



Moreover, using the freedom to choose the scale /i in eq. (3.10), the Z'th order 
coefficients (as far as they are known) can be brought down in magnitude below about 
(4tt)~ 1 , which means there is a fast decrease of terms in the perturbative series once 
< 1/3. This is shown in columns two and three of the figure. Unfortunately, the 
required scale \x is around a factor 4 or more below the mass of the quark. For the 
b-quark, a is rather large at that scale and the series is again unreliable. Only for even 
larger masses, say to* > 15 GeV, the asymptotic convergence of the series is noticeably 
better after adjusting the scale. More details are found in App. |A.2| Unfortunately we 
see no way out of the conclusion that for B-physics with a trustworthy error budget 
aiming at the few percent level, one needs a non-perturbative matching, even in the 
static approximation. 



We return to the full set of heavy- light flavor currents of Sect. 2.5 The bare fields 



satisfy the symmetry relations eq. (2.33). The same is then true for the RGI fields in 
static approximation. It follows that in static approximation the effective currents are 
given by 



Matching 35 



^jHQet = Cps (Mb/A _) zt% GI ( 9 o) A** , (3.36) 

7 fe HQET = C v (M h /A m ) Zf$ G1 (g ) Vt at , (3.37) 

V* QET = C PS (M h /A m ) Zf$ GI (g ) Z$?l(g ) V S ^ , (3.38) 

A* QET = C v (M h /A m ) Zf$ Gl (g )Z$f A (g ) Af at . (3.39) 

The factor RGi(ffo) is known as discussed in the previous lecture. Note that Z A r GI (<?o) 
is common to all (components of the) currents. Due to the HQET symmetries, there 
is one single anomalous dimension. A dependence on the different fields comes in only 
through matching, i.e. through the QCD matrix elements. In the above equations, chi- 



ral symmetry (of the continuum theory), eq. (2.30), has been used to relate conversion 
functions of axial and vector currents. 

Exercise 3.1 Pseudo-scalar and Scalar densities 
Start from the PCAC, PCVC relations in QCD 

<9 M (A R ) M = (fn b (p) + m 1 ( A t))P R ( M ) , (3.40) 
9m(Vr.)m = (m b (At) - m I (/x))5 R (/u) . (3.41) 

Replace all quantities by their RGI's. Take the matrix elements between vacuum and a 
suitable B-meson state to show that 

p hqet = _ Cps(Mb/A ) m* zttoMAT' , (3.42) 

Mb 

5 hqet = Cv(Mb/A _) !g Z$f A (g ) ZtnGi(go) V stat , (3.43) 
is valid up to terms of order 1/m. What happens if you choose a different matrix element? 

3.3.3 Applications 

As an application, we can now modify the scaling law for the decay constant to include 
renormalization and matching effects 



$rgi + 0(l/m) (3.44) 

(3.45) 



C PS (M h /A m ) 

/b _ y/m5 0p B (Mb/A m ) 



fu ^/m^Cp S (M c /A m )' 

where the latter equation is maybe stretching the applicability domain of HQET. 

Despite the discussion above, let us assume that the conversion functions C are 
known with reasonably small errors from perturbation theory. In this case, the knowl- 



edge of the leading term in expansions such as eq. (3.44) is very useful to constrain 
the large mass behavior of QCD observables, computed on the lattice with unphysical 
quark masses < m\,, typically mi, « m. c iiarm- (Such a calculation is done with a 
relativistic (Wilson, tmQCD, . . . ) formulation, extrapolating ami, — > at fixed m^.) 
As illustrated in Fig. |3.3[ one can then, with a reasonable smoothness assumption, 
interpolate to the physical point. 
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Fig. 3.3 Example of an interpolation between a static result and results with < mj,. The func- 
tion Cps is estimated at three-loop order. Continuum extrapolations are done before the interpolation 
I Blossicr et al., 2010cl. The point at 1/rom.pg = is given by 1"q &gx£l- This quenched computation 



is done for validating and demonstrating the applicability of HQET. 

Given the unclear precision of the perturbative predictions, the above interpolation 
method has to be taken with care. The inherent perturbative error remains to be 
estimated. 

The relation between the RGI fields and the bare fields has also been obtained for 
the two parity violating AB = 2 four fermion operators (Palombi et al, 200G Palombi 
et al, 20071) for N f = and N f = 2 (IDimopoulos et al, 20081). Their matrix elements, 



evaluated in twisted mass QCD will give the standard model B-parameter for B-B 
mixing. 

We now turn to the natural question whether one can directly compute the 1/m 
corrections in HQET, which will lead us again to the necessity of performing a non- 
perturbative matching between HQET and QCD. 



Exercise 3.2 Anomalous dimension 7 ma t c h 
Show that 



7match = -703* - [7l + 2&0Cl(l)]fl r * + . 



(3.46) 



where c\ is the 1-loop matching coefficient in the same scheme as 71. 



4 

Renormalization and matching at 
order 1/m 



4.1 Including l/m corrections 

We here work directly in lattice regularization. The continuum formulae are completely 
analogous. The expressions for 0ki n ,0 S pm & re discretized in a straight forward way, 



D k D k -> V^V fe , F kl -> F kl (4.1) 



with the latter given by the clover leaf representation, defined e.g. in (Liischer et al., 



1996). Of course other discretizations of these composite fields are possible. 

Apart from the terms in the classical Lagrangian, renormalization can in principle 
introduce new local fields compatible with the symmetries (but not necessarily the 
heavy quark symmetries which are broken by O sp i n , 0^ n ) and with dimension d op < 5. 
Also the field equations can be used to eliminate terms. With these rules one finds 
that no new terms are needed and it suffices to treat the coefficients of O sp i n , ©kin as 
free parameters which depend on the bare coupling of the theory and on m. 

The l/m Lagrangian then reads 

^h\x) = -Kin Okin(x) + CJ spin O spin (x)) . (4.2) 

Since these terms are composite fields of dimension five, the theory defined with a 
path integral weight (Jz?ii gn t collects all contributions of QCD with the heavy quark(s) 
dropped) 

W^nrqcd « exp(-a 4 ^[%t(i) + jS? h stat (z) + ^ l \x)]) (4.3) 

x 

is not renormalizable. In perturbation theory, new divergences will occur at each order 
in the loop expansion, which necessitate to introduce new counter-terms. The contin- 



uum limit of the lattice theory will then not exist(Thacker and Lepage, 19911. Since 
the effective theory is "only" supposed to reproduce the l/m expansion of the ob- 
servables order by order in l/m, we instead expand the weight W in l/m, counting 
w k i„ = 0(l/m) = o; spin , 

W^nrqcd -> WWt = exp(- a 4 ^[^ light (^) + J2f at (z)]) J 1 - a 4 ]T^ h (1) (*) 1 . 

X \ X ) 

This rule is part of the definition of HQET, just like the same step is part of Symanzik's 
effective theory discussed by Peter Weisz. 
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Let us remark here on the difference to chiral perturbation theory. In chiral per- 
turbation theory one computes the asymptotic expansion in powers of p 2 . Each term 
in the expansion requires a finite number of counter terms, since there are only a finite 
number of (pion) loops. The theory is thus renormalizable order by order in the ex- 
pansion. In NRQCD and HQET one expands in 1/m. At each order of the expansion 
an arbitrary number of loops remain, coming from the gluons and light quarks. In fact, 
we are even interested in more than an arbitrary number of loops: in non-perturbative 
results in a. 

NRQCD can then only be formulated with a cutoff and results depend on how the 
cutoff is introduced and on it's value. On the lattice, the cutoff is identified with the 
one present for the other fields, A cut ~ 1/a. Instead of taking a continuum limit, one 
then relies on physics results not depending on the lattice spacing (the cutoff) within 



a window ( Thacker and Lepage, 1991 1 



1/m < a < A QCD . [in NRQCD] (4.4) 

In HQET the discussion is rather simple, since the static theory is (believed to be) 
renormalizable; we will come to the renormalization of the insertion of shortly. 
Up to and including 0(l/m), expectation values in HQET are therefore defined as 

(O) = (0)stat + ^kinQ 4 y^(00kin(^))stat + 0J spin a 4 

slat 

x x 
= (O)stat + W kin (C) kin + W spin (0) S pi n , (4.5) 

where the path integral average 

(O)stat = \ f O exp(-a 4 yXSh&ix) + J2f at (z)]) (4.6) 

^ ./fields x 

is taken with respect to the lowest order action. The integral extends over all fields 
and the normalization Z is fixed by (l) s tat = 

In order to compute matrix elements or correlation functions in the effective the- 
ory, we also need the effective composite fields. At the classical level they can again 
be obtained from the Fouldy-Wouthuysen rotation. In the quantum theory one adds 
all local fields with the proper quantum numbers and dimensions. For example the 
effective axial current (time component) is given by 

^ ET (,) = Zl^ [Ar\ X ) + £ cf4\ X )\ , (4.7) 

i=l 

A£\x) = V^)l 757i (Vf- %Mx) , (4.8) 

A ( *\x) = -d i A? at , (4.9) 
where all derivatives are symmetric, 

1 A straight expansion gives e.g. o; kin a 4 5^2.(0 [Okin(^) — (Okin(is))stat)stat, but this just corre- 
sponds to an irrelevant shift of OkinM etc. by a constant. 
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$ = + Vf^c^ + V*), vf=i(Vi + vD, (4.10) 

and we recall ^ tat (ir) = - 0( a: )7i75'0h(^)- One arrives at these currents, writing down 
all dimension four operators with the right flavor structure and transformation under 
spatial lattice rotations and parity. The equations of motion of the light and static 
quarks are used to eliminate terms but heavy quark symmetries (spin and local flavor) 
can't be used since they are broken at order 1/mJ^] 

For completeness let us write down the other HQET currents: 

4 QET (*) = ^ QET ^ tat (*) + £^4%)] . (4-n) 

i=3 

A^(x) = ^)i 7fe757i (Vf- VfVh(x) , A£\x) = $(a:)i(Vg - Vl) 75 ^ h (x) , 

A { k 5 \x) = di (^(x) 7fc757 ^h(x)) , A[ 6 \x) = dkAt &t . 

The vector current components are just obtained by dropping 7 5 in these expressions 
and changing — > Cy^ . The classical values of the coefficients are Cj± — Cj^ = = 
— — while Cj^ = = 0. We note that with periodic boundary conditions in 
space we have 

a^^^^^^V^hW, a 3 £4 2) (x)=0, (4.12) 

XX X 

which for instance may be used in the determination of the B decay constant. 

Before entering into details of the renormalization, we show some examples how 
the 1/m-expansion works. 

4.2 1/m-expansion of correlation functions and matrix elements 

For now we assume that the coefficients 

0(1): 6m,Z^, 

(4.13) 

0(l/m) : w kin , w spin , c\> , 

are known as a function of the bare coupling g and the quark mass m. Their non- 
perturbative determination will be discussed later. 

The rules of the 1/m-expansion are illustrated on the example of C^ A ^(a;o), 



cq. (3.1 1. One uses eq. (4.5) and the HQET representation of the composite field 
eq. (4.7). Then the e xpecta tion value is expanded consistently in 1/m, counting pow- 
ers of 1/m as in eq. (4.13). At order 1/m, terms proportional to cjkin x Cp^ etc. are 
to be dropped. As a last step, we have to take the energy shift between HQET and 

2 An operator ^-Aff^ is included as a corresponding mass-dependence of Z^"^ ET . In practice, 
since <@C 1, and this term appears only at one- loop order, this dependence on the light quark 

m b 

mass can be neglected. 
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QCD into account. Therefore correlation functions with a time separation Xq obtain 
an extra factor exp(— Xom), where the scheme dependence of m is compensated by 
a corresponding one in Sm. Dropping all terms 0(l/m ) without further notice, one 
arrives at the expansion 



C*AA ( X o) 



(Z HQET )2 ^stat (xo) + Jl) C ff A ( XQ ) 
+ W k in C A a(x ) + W spin C A P A n (x ) 

mx °(ZI QET ) 2 Cft(x ) [l + 4 1} R$?{x ) 

+ ^kin RaK X o) + W spin R a P1 a( X o) 



(4.14) 



(4.15) 



with (remember the definitions in eq. (4.5)) 

Cl t A a i(x ) = a 3 ^« at (x)(4 1) (0))t> 8tat + a 3 ^2{A^(x)(Af at (0)y) 

X X 
X 

C A p A n (x ) = a 3 ^« at (^)« at (0))^ 



/spin • 



(2) 

The contribution of A vanishes due to eq. ( 4J.2 ). It is now a straight forward exercise 
to obtain the expansion of the B-meson mas^] 



m B = - lim d lnC? A D (x ) 

Xq— too 

= m bare - lim ab[lnCl t l t (x ) + c A 1) R$£{x ) + 



+ uj kin i? AA (x ) + w spin i? AA (a;o) ] 5m=0 



^kin = 



''"bare 

- lim 9 lnCI^Cxo) 



^stat 



kin 



pspin 



Sm=0 



- lim So^aaK), £ spm = - lim 9 i^ A n (x ) 



(4.16) 
(4.17) 

(4.18) 
(4.19) 

(4.20) 



Again we have made the dependence on Sm explicit through mbarc = m^+Sm and then 
quantities in the theory with Sm = appear. Note that the ratios i? AA (and therefore 
_E km , i<; s P m ) do not depend on 5m; the quantities _E km , E spin have mass dimension two 
and we have already anticipated eq. (4.25). 
The expansion for the decay constant is 

3 It follows from the simple form of the static propagator that there is no dependence on Sm except 
for the explicitly shown energy shift 8m. 



Renormalization beyond leading order 41 



/ B% Mi= lim {2exp(TO B ^o)C^ D (x )} 1/2 (4.21) 

Xq—^QO 

= Z* QET $ stat lim {l + \x [w kin £ kin + uj spin E s ^] 

a:o->oo 

+ ici 1) i?fl t (x ) + ^oW^^o) + H pin i? A p | n (z )} , (4.22) 
$ stat = lim {2cxp(^ stat a;o)C , i t l t (a;o)} 1/2 . 



Using the transfer matrix formalism (with normalization (B\B) = 2L 3 ), one further 
observes that (do it as an exercise) 

E kin = -^(B|a 3 ^0 kin (O,z)|i?) stat = -i(S|O kin (0)|fl) stat (4.23) 



E spin = -x(B\O spin (0)\B) stat , (4.24) 



0= lim doHfi?(xo). (4.25) 

Xq— >00 

As expected, only the parameters of the action are relevant in the expansion of hadron 
masses. 

A correct split of the terms in eq. (4.18 1 and eq. (4.22 ) into leading order and next 
to leading order pieces which are separately renormalized and which hence separately 
have a continuum limit requires more thought on the renormalization of the 1/m- 
expansion. We turn to this now. 

4.3 Renormalization beyond leading order 



For illustration we check the self consistency of eq. (4.14). The relevant question con- 
cerns renormalization: are the "free" parameters 5m . . . sufficient to absorb all 
divergences on the r.h.s.? We consider the term oc C^^(xo) since its renormalization 
displays all subtleties. As a first step we rewrite w k inCkin = (^"Or ^ n terms of a 
renormalized mass and the renormalized operator 

(0 kin ) R (z) = Z 0kin (O kin (z) + ^ if h (z)D i; h (z) + % fo,(*)Vh(*)) ■ ( 4 - 26 ) 

n a or 

The latter involves a subtraction of lower dimensional ones with dimensionless coeffi- 
cients Ci(go). The renormalization scheme for mg is irrelevant, as any change of scheme 
can be compensated by -Ze> kin , Cj whose finite parts need to be fixed by matching to 
QCD. We further expand 

(^HQETy = ( Z stat)2 + 2 Zf at Z^ /m) + 0(l/m 2 ) (4.27) 

which we will discuss more below. With these rules we then have 

(Zf at ) 2 u} kin Cfl(x ) = 5^a 7 G ( x > z ) + subtraction terms , (4.28) 



where 
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G(x,z) = ([Ar%(x) ([< at ] R (0))t (O kin ) R (z)\ . (4.29) 

\ / stat 

The subtraction terms are due to the lower dimensional operators with coefficients c\ 



and C2- Since we are interested in on-shell obscrvables (x > in eq. (4.14)), we may use 
the equation of motion DqiP^z) = to see that the ci-term does not contribut e, wh ile 



^J?/'h( z )' ! / , h(-2 : ), is equivalent to a mass shift. In the full correlation function eq. (4.14) it 
hence contributes to 6m which becomes quadratically divergent when the 1/m terms 
are included. 

While G(x, z) is a renormalized correlation function for all physical separations, its 
integral over z (or on the lattice the continuum limit of the sum over z) does not exist 
due to singularities at z — > and as z — > x. These contact term singularities can be 
analyzed by the operator product expansion. We discuss them first in the continuum 
and regulate the short distance region by just integrating for z 2 > r 2 with some small 
r. The operator product expansion then yields 



d 4 zG(x,z) (4.30) 

j 

[< a V(z) [di 1 - (A s tat (o)) f + d'UA^m^ + 4'(4 2) (o)) t 



stat 



up to terms which are finite as r — > 0. The coefficients d" in the operator product 
expansion have a further logarithmic dependence on rj^For (the continuum version 



of) eq. (4.29) we need r — > 0. In this limit short distance divergences emerge which 
have to be subtracted by counter-terms. In the lattice regularization, short distance 
singularities are regulated by the lattice spacing a and we have in full analogy 

[< at Mz) KEt^'lRiO)) 1 (^ in ) R (z)]) (4.31) 



[A^Mx) K- (^ tat (o)) t + 4(4 1) (o)) t + 4(4 2) (o)) t l 

CL I stat 

up to terms which have a continuum limit a — > and up to the singular terms originat- 
ing from z R3 x . The coefficients di contain a logarithmic dependence on a. Treating 
the singular terms at z ~ x in the same way and noting that the term with A Q (0) 
vanishes upon summation over x we find 

Zl tat [ d x - Cff (x ) + d 2 Cff (xo) } (4.32) 
a 



for the contact term singularities in eq. (4.28 1. These are absorbed in eq. (4.14 ) through 



counter-terms contained in Z^® ET and c^, 

97 (1/m) - dl I r (1) - d ' 2 I (dW) 

The change from d\ to di is due to the use of the equation of motion above. This step 
is valid only up to contact terms, resulting in the shift d' — s- d. The ellipses contain 
the physical, finite 1/m terms. 



4 We have written down the integrated version, since then a smaller number of operators can 
appear and we are ultimately interested in the integral. 



The need for non-perturbative conversion functions 43 



We now comment further on the expan sion eq. (4.27|. Our discussion shows that 



the quadratic term (Z^" 1 ^) 2 in eq. (|4.27| must be dropped] otherwise an uncanceled 



l/(a 2 m 2 ) divergence remains. As we have seen there is no l/(a 2 m 2 ) in CY^(xo) and 



the other pieces in eq. (4.141 are less singular. This is just a manifestation of the 



general rule of an effective field theory that all quantities are to be expanded in 1/m 
whether they are divergent or not. With this rule the various HQET parameters can 
be determined such that they absorb all divergences. 

The lesson of our discussion is that counter-terms with the correct structure are 
automatically present because in the effective theory all the relevant local composite 
fields are included with free coefficients. These free parameters may thus be chosen such 
that the continuum limit of the HQET correlation functions exists. Finally, their finite 
parts are to be determined such that the effective theory yields the 1/m expansion of 
the QCD observables. 



4.4 The need for non-perturbative conversion functions 

An important step remains to be explained: the determination of the HQET param- 
eters. As discussed in Sect. [3] at the leading order in 1/m, this can be done with the 
help of perturbation theory for conversion functions such as Cps- However, as soon 
as a 1/m correction is to be included, the leading order conversion functions have 
to be known non-perturbatively. This general feature in the determination of power 



corrections in QCD is seen in the following way. Consider the error made in eq. (3.9), 
when the anomalous dimension has been computed at / loops and C ma t c h at I 
order. The conversion function 



1 loop 



C PS = exp 



9* "f X^ 

dx 



match ™2/ 



p{x) 



+ A(Cps) 



(4.34) 



is then known up to a relative error 



A(C; 



PS) 



Cps 



oc [g 2 {m) 



2b a ln(m/A QC D) 



A 



QCD 



(4.35) 



As m is made large, this perturbative error becomes dominant over the power correc- 
tion one wants to determine. Taking a perturbative conversion function and adding 
power corrections to the leading order effective theory is thus a phenomenological ap- 
proach, where one assumes that for example at the b-quark mass, the coefficient of the 
[<? 2 (m,b)] z term (as well as higher order ones) is small, such that the A/m^ corrections 

5 It is convenient to avoid the multiplication of 1/m terms explicitly by a choice of observables, 
for example 

4 = ln(/ Bv ^B) = m(z£ QET ) + ln(<t> stat ) + lim { ±x uj kin E kin + iw Mn R^l(xo) + ...}, 



7 (l/m) 



) = ln(Zi tat ) + 



ln(Z A 

In this convention all 1/m-terms appear linearly. 



'Zstat 
Z A 



ln(3S at ) + lHZA)] 1/m 
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dominate. In such a phenomenological determination of a power correction, its size 
depends on the order of perturbation theory considered. A theoretically consistent 
evaluation of power corrections requires a fully non-perturbative formulation of the 
theory including a non-perturbative matching to QCD. Note that the essential point 



of Eq. (4.35 ) is not the expected factorial growth of the coefficients of the perturbative 
expansion. Rather it is due to the truncation of perturbation theory as such. Of course 
a renormalon-like growth of the coefficients does not help. 

The foregoing discussion is completely generic, applying to any regularization. 
When we define the theory on the lattice, there are in addition power divergences, 



e.g. in eq. (4.31). It is well known that they have to be subtracted non-perturbatively 



if one wants the continuum limit to exist. 



4.5 Splitting leading order (LO) and next to leading order (NLO) 

We just learned that the very definition of a NLO correction to /b means to take 
eq. (4.22 1 with all coefficients Z^ ET . . . determined non-perturbatively. We want 



to briefly explain that, as a consequence, the split between LO and NLO is not unique. 
This is fully analogous to the case of standard perturbation theory in a, where the 
split between different orders depends on the renormalization scheme used, and on the 
experimental observable used to determine a in the first place. 



Consider the lowest order. The only coefficient needed in eq. ( 4.22 1 is then Z^ ET = 
Cps-^A rgp K has to be fixed by matching some matrix element of Ag tat to the matrix 

element of Aq in QCD. For example one may choose (S'|Ao|0), with \B') denoting some 
other state such as an excited pseudo-scalar state. Or one may take a finite volume 
matrix element defined through the Schrodinger functional as we will do later. Since the 
matching involves the QCD matrix element, there are higher order in 1/m "pieces" in 
these equations. There is no reason for them to be independent of the particular matrix 
element. So from matching condition to matching condition, Cps^a rgi determined 
at the leading order in 1/m differs by 0(AQCD/wb) terms. 

The matrix element /b in static approximation inherits this 0(AqcE)/mb) ambigu- 
ity. These corrections are hence not unique. Fixing a matching condition, the leading 
order /b as well as the one including the corrections can be computed and have a con- 
tinuum limit. Their difference can be defined as the 1/m correction. However, what 
matters is not the ambiguous NLO term, but the fact that the uncertainty is reduced 
from 0(AQCD/wb) in the LO term to 0(Aq CD /m^) in the sum. 

The following table illustrates the point explicitly. 

Observables (^|^|0) (B'|^J|0) {B"\A\\Q) 

matching condition * 

error in HQET result 0(A/m b ) 0(A/m b ) 



matching condition 

error in HQET result 0(A/m b ) 0(A/m b ) 

matching condition * 

error in HQET result 0(A/m b ) 0(A/m b ) 
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As a consequence, there is no strict meaning to the statement "the 1/m correction to 
/ B is 10%". 

4.6 Mass formulae 

Often cited mass formulae are 

m| v ee ][m B +3m B -] = m b + A+-— A x +0(l/m?) (4.36) 

4 Z?7lb 

2 o 

Am B = m B * - m B = A 2 + 0(l/mg) (4.37) 

m b 

with (ignoring renormalization) 

Ai = (B\O kin \B) , A 2 = l(B\O spin \B) . (4.38) 

The quantity A is termed "static binding energy" . Also here, depending on how one 
formulates the matching condition which determines mb, one changes A by a term of 
order Aqcd, e -g- one may define A = 0. Similarly, the kinetic term Ai/(2mb) has a 
non-perturbative matching scheme dependence of order Aqcd and thus Ai itself has a 
matching scheme dependence of order mb- The situation for A is similar to the gluon 
"condensate" . The non-perturbative scheme dependence has the same size as the gluon 
"condensate" itself. In contrast, A 2 is the leading term in the 1/m expansion and does 



not have such an ambiguity. We refer also to the more detailed discussion in ( Sommer, 



2006) 



4.7 Non-perturbative determination of HQET parameters 

We close our theoretical discussion of HQET by stating the correct procedure to de- 
termine the Ahqet parameters in the effective theory at a certain order in 1 /m. One 
requires 

^ QCD (m) = <i>f QET (m,a), i = 1 . . . iV HQ ET , (4.39) 

where the m-dependence on the r.h.s. is entirely inside the HQET parameters. On the 
l.h.s. the continuum limit in QCD is assumed to have been taken, but the r.h.s. refers 
to a given lattice spacing where it defines the bare parameters of the theory at that 
value of a. We emphasize that as this matching has to be invoked by numerical data, it 
is done at a given finite value of 1/m. Carrying it out with just the static parameters 
defines the static approximation etc. 



As simple as it is written down, it is non-trivial to implement eq. (4.39) in practice 
such that 

1) the HQET expansion is accurate and one may thus truncate at a given order, 

2) the numerical precision is sufficient, 

3) lattice spacings are available for which large volume computations of physical 
matrix elements can be performed. 

In the following section we explain how these criteria can be satisfied using Schrodinger 
functional correlation functions and a step scaling method. The first part will be a test 
of HQET on some selected correlation functions. This establishes how 1) and 2) can 
be met. We can then explain the complete strategy which also achieves 3). 
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4.8 Relation to RGI matrix elements and conversion functions 



The matching equations eq. (4.39) provide a definition of all HQET parameters, in 



principle at any given order in the expansion. If considered at the static order, it also 
provides the renormalization of the static axial current, which we discussed at length 
m Sect. [3} The relat ion between the two ways of parametrizing the current in static 
approximation are 



7 HQET 



0(l/m,a) 



While eq. (4.39) is a matching equation determining directly the product Z' 



(4.40) 

A,RGI C 'PS> 



the r.h.s. separates the problem into a pure HQET problem, the determination of the 
RGI operator, and a pure QCD problem, the the "anomalous dimension" 7 ma t c h, see 
eq. (3.201. Note that in this simple form, such a separation is only possible at the 



lowest order in 1/m. 

Since the breaking of spin symmetry is due to a single operator at order 1/m, there 
is also an analogous representation of w sp in- We refer the interested reader to (Guazzini 



et a/., 2007). 



5 

Non-perturbative HQET 



After our long discussion of the theoretical issues in the renormalization of HQET, we 
turn to a complete strategy for the non-perturbative implementation. To this end the 
three criteria in Sect. 4.7 have to be fulfilled. Establishing 1) is equivalent to testing 



HQET. We therefore start with such a test. Item 2) has to do with finding matching 
conditions sensitive to the 1/m-suppressed contributions. For this purpose we then 
expand a little on correlation functions in the Schrodinger functional before coming to 
a full description of the matching strategy. 



5.1 Non-perturbative tests of HQET 

Although it is generally accepted that HQET is an effective theory of QCD, tests of this 
equivalence are rare and mostly based on phenomenological analysis of experimental 
results. A pure theory test can be performed if QCD including a heavy enough quark 
can be simulated on the lattice at lattice spacings which are small enough to be able 



to take the continuum limit. This has been done in the last few years (Heitger et at, 



2004 Delia Mortc et at, 20081 and will be summarized below. 



We start with the QCD side of such a test. Lattice spacings such that am^ <C 1 can 
be reached if one puts the theory in a finite volume, L 3 x T with L, T not too large. 
We shall use T = L. For various practical reasons, Schrodinger functional boundary 
conditions are chosen. Equivalent boundary conditions are imposed in the effective 
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Fig. 5.1 Testing eq. thr ough numerical simulations in the quenched approximation and for 

L pa 0.2 fm I Heitger et al, 2004). The graph uses notation V^ CD = Kps). The physical mass of the 
b-quark corresponds to z ss 6. Two different orders of perturbation theory for Cps are shown. 
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theory. As in Sect. 3.3 we consider the ratio Y® CB (8, rriR, L) built from the correlation 
functions /a and f\. 



It can be written as 

v QCD, fl n (Q(L)\A \B(L)) 

|B(L)) = e- iH / 2 |^ B (i)) , |fi(L)> = c-™/ 2 ML)) , 

in terms of the boundary states \<Pb(L)) , \(po(L)). Expanded in energy eigenstates 
with energies E n > me in the B-sector and energies E n in the vacuum sector, we have 

\B(L)) = e' LE " /2 (n, B\^(L)) \n, B) (5.2) 

n 

£ e- LE "/ 2 (n,B\ l p B (L))\n 1 B)+0(e- k / 2 ), (5.3) 

n | E n — m-B<k/ L 

|n(£)) = e- LS " /2 (n, 0\ML)) |«, 0) , (5.4) 
~ £ e- i ^/ 2 <n,0|^ (i)>|n,0) + O(e- fc / 2 ), (5.5) 

n | E n <k/L 

which shows that only energy eigenstates with E n — Eq = 0(1/ L) contribute signifi- 
cantly. For z = LMb 3> 1, HQET will thus describe the correlation functions and the 
ratio We come to the conclusion that 

Y^ CD (8,m R ,L) = C PS (M b /A)X RGI + 0(1/ z) , z = M h L, (5.6) 
X RGI = lim Zt% GI (9o)^^ (5.7) 

and similarly for other observables. Note that one could also just argue that the only 
relevant scales are L, A, m^. Therefore with L ss 1/A there is a A/mb ~ 1/z expansion. 



Of course relations such as eq. (5.6 1 are expected after the continuum limit of both 
sides has been taken separately. For the case of Y® CB , this is done by the following 
steps: 

• Fix a value uq for the renormalized coupling g 2 (L) (in the Schrodinger functional 



scheme) at vanishing quark mass. In (Heitger et al., 2004) uq was chosen such 
that L w 0.2 fm. 

For a given resolution L/a, determine the bare coupling from the condition g 2 (L) = 



Uq. This step is well known by now (Capitani et al., 1999). 

Fix the bare quark mass m q of the heavy quark such that LM = z using the known 
renormalization factors Zyi, Z in M = Zy^Z (1 + ab m m q ) m q , where Z, Zyi, b m a re 
all known non-perturbatively ( Guagnelli et al., 2001||Della Morte et al., 2007a ). 



• Evaluate Y® CD and repeat for better resolution a/L. 

• Extrapolate to the continuum as shown in Fig. |5.1( left. 
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-1.284 - 
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Fig. 5.2 Continuum extrapolation of Xrgi (Heitger et al, 2004). 

In the effective theory the same steps are followed. As a simplification, no quark 
mass needs to be fixed and the continuum extrapolation is much easier as illustrated 
in Fig. [IT2] 

The comparison of the static result and the relativistic theory, Fig. |5.1| looks rather 
convincing^ but we note that the b-quark mass point is 1/z = 1/zb ~ 0.17, where 
1/z 2 terms are not completely negligible. The displayed fit has a 8% contribution by 
the 1/z term and a 2% 1/z 2 piece. 

For a precision application (Sect. 5.3) it is thus safer to have L>0.4fm instead of 
the L = 0.2 fm chosen in the first test, reducing 1/z 2 by a factor four. For L ss 0.5 fm 



we show two different examples, Fig. 5.3 Fig. 5.4 which involve 



ki(0) 



6L 6 



£ (C,'(u) 7 *Cb(vRb(y)7 fc Ci(z)) 

u,v,y ,z,fe 



(5.8) 



in addition to the previously introduced correlation functions. The considered combi- 
nations are 



Ri 
Ri 



i 1,1 /i(fli)fci(3i) 3 

4 1 l Jl(*2)*l(fc) 3 



4 Ui 



(5.9) 
(5.10) 



Their HQET expansion contains no conversion functions at leading order and they are 
thus free of the associated perturbative uncertainty. While R\ has a finite static limit, 
R\ vanishes as z — ¥ oo due to the spin symmetry. The expected HQET behavior is 
confirmed with surprisingly small 1/z 2 corrections for a charm quark. The quadratic 
fits in 1/z displayed in the figures are not constrained to pass through the separately 
displayed static limit. 



1 Note that the comparison Fig. 5.1 has to be taken with a grain of salt due to the perturbative 



uncertainty in Cpg discussed in Sect. |3..'Oj 
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Fig. 5.3 The logarithmic ratio Ri for different pairs (61,62) with Nf = 2 flavors and for 
L = T pb 0.5 fm (Delia Morte et al, 2008[ > with N f = 2. The value of 1/z for charm and bottom 
quarks are indicated by the vertical bands. 
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Fig. 5.4 The logarithmic ratio R\ for different values 80 with N{ = 2 flavors and for L = T £3 0.5 fm 
(Delia Morte et al, 2008] ) with Nf = 2. 

5.2 HQET expansion of Schrodinger functional correlation 
functions 



In complete analogy to the case of a manifold without boundary we can write down the 
expansions of the Schrodinger functional correlation functions to first order in 1/m: 
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[/ A ] R = z^ ET z Ch z c e- m ^« {/i tat + c^/ir + cwJr+<*wf n } #.n) 

[/i] R = Z 2 Ch Z 2 c e- m — T {/f at + w kin A kin + ^ P i„A spin } , (5.12) 
[fci] a = ZlZ 2 c e~ m — T {/f at + Wki „A kin - ^spin/i pin } ■ (5.13) 



Apart from 



^(«o,e) = -y E (4 1) WCh(y)7 5 C 1 (z)) (5.14) 



the labeling of the different terms follows directly the one introduced in eq. (4.5). 
The relation between the 1/m terms in /i and k\ is a simple consequence of the spin 
symmetry of the static action, valid at any lattice spacing. A further simplicity is that 
no 1/m boundary corrections are present. Potential such terms have dimension four. 
After using the equations of motion, only one candidate remains, which however does 
not contribute to any correlation function^] 

5.3 Strategy for non-perturbative matching 

After the tests of HQET described above, it is clear how one can non-perturbatively 
match HQET to QCD. Consider the action as well as Ao (just at p = 0) and denote the 
free parameters of the effective theory by u>i , i = 1 . . . Ahqet- In static approximation 
we then have 

w stat = ( m stat e ; [m(ZA) ]stat )t ^ A/jjQET = 2 (5.15) 

and including the first order terms in 1/m together with the static ones, the HQET 
parameters are 



HQET = ( mbarc ^ ln(Z HQE T) _ ^1) > ^ ^ ^ y = g _ (g lg) 



The pure 1/m parameters may be defined as a/ 1 /™) — W H Q ET — w stat j with all of them, 

e.g. also m^™ 1 ', non-zero. In fact our discussion of renormalization of the 1/m terms 

shows that m^.™' 1 diverges as l/(a 2 m). 
With suitable observables 

$,(ii,M,o), i = l...JV H QBT, 

in a Schrodinger functional with L = T = L\ ~ 0.5 fm, we then require matching 

$ i (£ 1 ,M,a) = $9 CD (£ 1 ,M,0), i = l...JV H QET. (5.17) 



2 In the notation of (Liischer et aL, 19961 it reads ph(x)7fc.D fcPh( x ) a * = 0. Such a term does 
not contribute to any correlation function due to the form of the static propagator. 
^Recall that observables without a superscript refer to HQET. 
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Note that the continuum limit is taken in QCD, while in HQET we want to extract 
the bare parameters of the theory from the matching equation and thus have a finite 
value of a. It is convenient to pick observables with HQET expansions linear in taj, 

$(L, M, a) = T](L, a) + <f>{L, a) oj(M, a) , (5.18) 

in terms of a -/Vhqet x -/Vhqet coefficient matrix <j>. A natural choice for the first two 
observables is 

$1 = LT P = -Ld H-f A {x Q )) X0=L /2 L ~°° Lm B (5.19) 

$ 2 = ln(Z A ^) ^X^/bv^, (5.20) 
v h 

since in static approximation these determine directly u>i and w 2 . We will introduce 
the other <3>i later. The explicit form of 77, (f> is 

/r stat \ / L . . A 

T, = Ca , 0-01... (5.21) 

with 

j*stat 

T s ^ = -Ld ln(f^(x )) X0=L/2 , U=H-f=)- (5-22) 

V h 

In static approximation, the structure of the matrix cj) is perfect: one observable de- 
termines one parameter. This is possible since there is no (non-trivial) mixing at that 
order. 

Having specified the matching conditions, the HQET parameters 0Ji{M,a) can be 
obtained from eqs.( 5.17|5.18 1, but only for rather small lattice spacings since a reason- 



able suppression of lattice artifacts requires L\/a = O(10) and thus a = 0(0.05 fm). 
Larger lattice spacings as needed in large volume, can be reached by adding a step 



scaling strategy, illustrated in Fig. |5.5| Let us now go through the various steps of this 
strategy. 

(1) Take the continuum limit 

$9 CD (£ 1 ,M,0)= lim $^ CD (Li, Af, a) . (5.23) 

a/Li— >0 

This is similar to the HQET tests and as we saw there, it requires Li/a = 20 . . . 40 , 
or a = 0.025 fm. . .0.012 fm. 



(2a) Set the HQET observables equal to the QCD ones, eq. (5.17) and extract the 
parameters 

w(M, a) = <tr\Li,a) [$(L lf M, 0) - r,(L u a)] (5.24) 

^7 1 $i(L 1 ,M,0)-r stat (Li,a)\ 

$ 2 (£i,AT,0)-CA(£i,a) ■ (5.25) 

The only restriction here is L%/a 1, so one can use L\ja — 10 ... 20 , which 
means a = 0.05 fm . . . 0.025 fm. 
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HQET 




Fig. 5.5 Strategy for non-perturbative HQET | |Blossier et ah, 20 10b|. Note that in the realistic 
implementation| |Blossier et al., 20106^ finer resolutions are used. 



(2b.) Insert Q into $(£ 2 ,M, a) 



$(L 2 ,M,0) = lim {r/(L 2 ,a) + 0(L 2 ,a)w(M,a)} 

a/L 2 — i-0 



(5.26) 



/L 2 r stat (L 2 , o) + (Li , M, 0) - L 2 r stat (L u a) 
lim CA(i2,a) + $ 2 (ii,M,0)-CA(£i,a) 

%[r stat (L 2 ,a) -r stat (L!,a)]\ M, 0) N 

lim I CA(£2,a)-C A (£i,a) + $ 2 (£i,M,0) 

a/ L 2 ->0 



finite HQET SSF's 



QCD, mass dependence 



In the last line we have identified pieces which are separately finite. This step 
can be done as long as the lattice spacing is common to the n 2 = L 2 /a and 
ni = Li/a-lattices and 



s = L 2 /Li = n 2 /n 1 



(5.27) 
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is kept at a fixed, small, ratioj^] 
(3.) Repeat (2a.) for L\ — > L 2 : 

w(M, a) = (p- 1 ^, a) [$(L 2 , M, 0) - T)(L 2 , a)] . (5.28) 

With the same resolutions -L 2 /a = 10 ... 20 one has now reached a = 0.1 fm . . . 0.05 fm. 
(4.) Finally insert 10 into the expansion of large volume observables, e.g. 

m B = wi + £ stat . (5.29) 

In the chosen example the result is the relation between the RGI b-quark mass and 
the B-meson mass me- It is illustrative to put the different steps into one equation, 



(5.30) 



lim[£; stat -r stat (L 2 ,a)l a =0.1fm . . . 0.05fm [S 4 ,S 5 ] 

a— !-0 

+ lim[r stat (L 2 ,a) -r stat (L!,a)] a =0.05fm . . . 0.025fm [S 2 ,S 3 ] 
+ — lim M b , a) a =0.025 fm . . . 0.012 fm [Si] . 

_Ll a— >0 

We have indicated the lattices drawn in Fig. |5.5| and the typical lattice spacings of 
these lattices. The explicit expression for the decay constant in static approximation 
is even more simple; write it down as an exercise! 

So far we have spelled out only those observables which are needed in the static 
approximation. The following heuristics helps to find observables suitable for the de- 
termination of the 1/m-terms. Recall that 9^0 means |(Vj + V*) ~ i9/L (acting 
onto a quark field) when the gauge fields are weak, as is the case in small volume. 
Hence, expanding in 1/m 



4> 3 (L, M. a) = ^1 ~ . . . + c£> [0 2 - (5.31) 



for weakly coupled quarks. In the same way the combination (recall eq. (5.9|) 

$ 4 (L, M, a)=R 1 = i?f at + u Mn Rf n 

has a sensitivity to Wkin of i?i' n oc Q\ — Q\ while in the specific linear combination of 
/i and k\ which form R\ the parameter cj sp i n drops out. Finally the choice 

$ 5 (L, M, a) = R 1 = tJ spia Rf n (5.32) 

allows for a direct determination of w sp i n . These choices leave relatively many zeros in 
the matrix which has a block structure, 

C B\ (C- x -C- X BA- X \ „ fLO 



0AJ> r A ~' J' ^^V 01 /' (5 ' 33) 

The listed observables <I>i have been shown to work in practice, i.e. in a numerical 
application (Blossier et al., 20106K 

4 A fixed ratio s ensures that the cutoff effects are a smooth function of a/Li. 
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5.4 Numerical computations in the effective theory 

Before showing some results, we should briefly mention that it is not entirely straight 
forward to obtain precise numerical results in the effective theory. The reason is a 
generically rather strong growth of statistical errors as a function of the Euclidean 
time separation of the correlation functions. Two ideas help to overcome this problem. 
We sketch them here; more details are available in the cited literature. 

5.4.1 The static action 



Consider a typical two-point function, for example eq. (3.5 1. At large time it decays 



exponentially and so does the variance. Setting 5m — the decay of the signal is 

C(x ) ~ e - E ^ tXo , (5.34) 

while the variance decays with an exponential rate given by the pion mass. Thus the 
noise-to-signal ratio for the B-meson correlation function behaves as 

i? NS OC e [^t-m r /2]x (5 35) 

The self energy of a static quark is power divergent, in particular in perturbation 
theory 



Q r (i) +0 (a°)) .g 2 + O(g 4 ). (5.36) 



This divergence yields the leading behavior of eq. (5.351 for small a. It is potentially 
dangerous since we are interested in the continuum limit. The scale of the problem 
can be reduced considerably by the replacement 

I7(a;,0) -> WHYPifoO), (5.37) 

in the covariant derivative Vo in the static action. Here WuYPi is a so-called HYP- 
smeared link. Table [5~T] shows how the self energy is reduced for two choices of W^hypi- 



Sf H 0.16845(2) 0.68(9) 
^ YP1 0.04844(1) 0.44(2) 
S* YP2 0.03523(1) 0.41(1) 

Table 5.1 One loop coefficients r' 1 ', eq. |5,36[ l and non-perturbativc values for ai? s tat at 
/3 = 6/gg = 6 and a (quenched) light quark with the mass of the strange quark. "EH" refers to 
Eichten-Hill, i.e. W(x,0) = U(x,0), while "HYP1,HYP2" are two versions of HYP-smearing ( |Hasen-| 
|fratz and Knechtli, 20uT)|Della Morte et al, 2005| l 



It is mandatory to check that such a change of action does not introduce large 
cutoff effects. This was done for single smearing in (Delia Morte et at, 2005): the 



points with smallest error bars in Fig. |2.1| are for these actions. We expect that large 
cutoff effects would however appear if smearing was repeated several times. 
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Fig. 5.6 Numerical solution of the equation for | |Blossier et al., 20106[ l made dimensionless by 
multiplication with L2. The figure use 
<E>2 in the figure is $1 in our notation. 

5.4.2 Generalized Eigenvalue Method 



For the numerical evaluation of matrix elements such as q> stllt : eq. (1.57), or of energy 



levels it is advisable to use an improvement over the straight forward formula eq. (3.7) 



The reason is as follows. Let us label the energies in the sector contributing to a given 
correlation function by E n , n = 1,2,3. Then there are corrections to the desired 
ground state matrix element due to excited state contaminations of order e" x ° A and 
A = E 2 — Ei . From an investigation of the spectrum in the B- meson sector one finds 
numerically A s» 600 MeV and thus A xq « Sx$ /fm. The suppression of excited state 
contaminations is then not necessarily small enough for xq ~ lfm but using eq. (3.7) 
beyond xq ~ lfm is very difficult because statistical errors grow quite rapidly with xq. 
A considerable improvement is achieved if one considers the generalized eigenvalue 



problem (GEVP)( Michael and Teasdale, 1983 Luscher and Wolff, 1990 Blossier et al, 



2009). It uses additional information in the form of a matrix correlation function 



formed from N different interpolating fields on one time slice and the same interpolat- 
ing fields on another time slice. When this matrix correlation function is analyzed in a 



specific way, described in (Blossier et al., 2009), one can prove that a much larger gap, 
A = En+i — Ei appears for the dominating correction terms due to excited states. 
These then disappear much more quickly with growing time. 

The GEVP is straight forwardly applicable to HQET, order by order in 1/m. The 
precision of the numerical results that we show below is largely due to this method, 
together with the use of HYP1/2 actions. 



5.5 Examples of results 



We now discuss a few numerical results (Blossier et al., 2010& Blossier et al, 2010a 



Blossier et al., 2010c) in order to give an indication of what can be done at present. 
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LO (static) 


NLO (static + 0(l/m)) 






(0i,0a) = (0,0.5) 


(0i,0 2 ) = (0.5,1) ( 


01,02) = (0,1) 


#0 = 


17.1 ±0.2 


17.1 ±0.2 


17.1 ± 0.2 


17.1 ±0.2 


O = 0.5 


17.2 ±0.2 


17.2 ±0.2 


17.2 ±0.2 


17.1 ±0.2 


0o = l 


17.2 ±0.2 


17.3 ±0.3 


17.3 ±0.3 


17.3 ±0.3 



Table 5.2 Dimensionless b-quark mass, rgMb, obtained from the B B meson mass, for differ- 



ent values of 8i . 
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a 2 [fin 2 ] a 2 [fm 2 ] 

Fig. 5.7 Continuum extrapolations in HQET. Left: 5> H Q ET = / Bs ^r==/C P s (diamonds) in HQET 
with 1/m corrections includedl |Blossier et al., 2010c^ and its static limit "1?rgi (circles). The value of 
Cps does not depend on the lattice spacing. It renders the two quantities directly comparable. Right: 
pseudo scalar energy levels l |Blossier et al., 2010a[ |. From bottom to top: 2s — Is splitting static, 2s - 
Is splitting static + 1/m, 3s - Is splitting static. 
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Fig. 5.8 Static results together with results with < mj, and an HQET computation with 1/m 
corrections included. Continuum extrapolations are done before the interpolation ( |Blossier et al.,] 
2010 c|' Cps i s evaluated with the three-loop approximation of 7match- 
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The graphs and numbers are for the quenched approximation (the light quark is a 
strange quark) but these computations are also on the way for dynamical fermions. 
The statistics employed in the quenched approximation is rather modest: only 100 
configurations were analyzed. One can easily use a larger number, even with dynamical 
fermions. We skip numerical details in the following discussion. 

As a first step, one wants to fix the b-quark mass. This is done through eq. (5.30) 
and its 1/m corrections. Its graphical solution is illustrated in Fig. 5.6 where all plotted 
numbers originate from prior continuum extrapolations. The resulting mass of the b- 
quark is displayed in Table |5.2| Observe that it depends very little on the matching 
condition, i.e. the choice of 80,61,62 and moreover the 1/m corrections are small. 

Next we look at the lattice spacing dependence of the decay constant. For the 



results including 1/m corrections no significant dependence on a is seen in Fig. 5.7 



despite a good precision of about 2%. In static approximation, discretization errors 
are visible but small. Table [53| lists the B s decay constant using tq = 0.5fm to convert 
to MeV for illustration. The actual number is affected by an unknown "quenching 
effect" and thus not so important. It is more relevant to observe the precision that can 
be reached with just 100 configurations and how the spread in the numbers in static 
approximation is reduced when the 1/m corrections are included. 





LO (static) 


NLO (static + 0(l/m)) 






(0i, 2 ) = (0,0.5) 


(0i,0a) = (0.5,1) ( 


0!,0 2 ) = (O,1) 


0o 


= 


233 ±6 


220 ±9 


218 ±9 


218 ±9 


6o-^ 


= 0.5 


229 ±7 


221 ±9 


219 ±8 


219 ±9 


9o 


= 1 


219 ±6 


223 ±9 


221 ±8 


222 ±8 



Table 5.3 Pseudo-scalar heavy-light decay constant /b b in MeV, for different values of 6i. 



Further, the comparison with results in the charm mass region, Fig. 5.8 



to indicate that the 1/m expansion works very well even for charm quarks. This is 
a bit surprising and certainly requires further confirmation. Note also that this com- 
parison makes use of the perturbatively evaluated Cps whose intrinsic uncertainty 
due to perturbation theory is difficult to evaluate. Of course this uncertainty does 
neither affect the non-perturbatively computed static value at 1/ (romps) = 0, nor 
/b s v /m B s computed with 1/m corrections at the mass of the b-quark, corresponding 
to l/(r m PS ) R3 0.07. It only affects the comparison to the results for l/(r o mps)^0.15 
since only for the purpose of this comparison the logarithmic mass dependence de- 
scribed by Cps has to be divided out. 

Finally we show some results concerning the spectrum. The splitting between radial 
excitations in the pseudo-scalar sector is displayed in the right part of Fig. |5.7| As 
throughout in our results, the 1/m-correction is rather small. 



5.6 Perspectives 

Meanwhile it has been established that HQET with non-perturbatively determined 
parameters is a precision tool. However, we are still at the beginning concerning appli- 
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cations. Results for the quantities shown here will be available for N { = 2 dynamical 
fermions rather soon. But there are many more applications which remain unexplored 
and open interesting avenues of research for the future. 



Appendix A 



A.l Notation 

A. 1.1 Index conventions 

Lorentz indices fj,, v, . . . are taken from the middle of the Greek alphabet and run from 
to 3. Latin indices k,l, . . . run from 1 to 3 and are used to label the components 
of spatial vectors. For the Dirac indices capital letters A, B, . . . from the beginning 
of the alphabet are taken. They run from 1 to 4. Color vectors in the fundamental 
representation of SU(7V) carry indices a, (3, . . . ranging from 1 to N, while for vectors 
in the adjoint representation, Latin indices a,b, . . . running from 1 to TV 2 — 1 are 
employed. 

Repeated indices are always summed over unless otherwise stated and scalar prod- 
ucts are taken with Euclidean metric. 

A. 1.2 Dirac matrices 

In the chiral representation for the Dirac matrices, we have 




(A.l) 



The 2x2 matrices e u are taken to be 



= -1 



e k = -ia k , 



(A.2) 



with a k the Pauli matrices. It is then easy to check that 

7/^ = 7^ {7^,74 = 2( W 
Furthermore, if we define 75 = 70717273, we have 



(A.3) 




(A.4) 



In particular, 75 = 7J and 75 = 1. The hermitian matrices 
are explicitly given by (<Ji<Jj — itijkOk) 



(A.5) 




(A.6) 



where e^-fc is the totally anti-symmetric tensor with e 12 3 = 1. 
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In the Dime representation we have 



- - (,: -7) . » - c -1) ■ <-» 

75 = ( ? n ) 5 CT y = f °. ) = ^ ( A - 8 ) 



A. 1.3 Lattice conventions 

Ordinary forward and backward lattice derivatives act on color singlet functions f(x) 
and are defined through 

d ll f(x) = ±[f(x + afi)-f(xj\, 

dlf{x)= l a [f{x)~ f(x-ajij\, (A.9) 
where // denotes the unit vector in direction fj,. We also use the symmetric derivative 

^ = i(3„ + d*). (A.10) 
The gauge covariant derivative operators, acting on a quark field ip(x), are given by 

V^{x) = - [\„U(x, fi)ip(x + afi) - 1>{xj\ , (A.ll) 

V*^{x) = - a [tl>(x) - X-'Uix - aft, - oA)] , (A.12) 

with the constant phase factors 

A M = e ia6 ^ L , 9 = 0, -ir<9 k <ir, (A.13) 



explained in Sect. 2.6 The left action of the lattice derivative operators is defined by 
^(x)V M = - [^(x + aftUix,^)- 1 ^ 1 -i>(x)] , (A.14) 

ip(x)X7* = - \ip(x) - ip(x - afi)U(x - afi,fJ,)Xu] • (A. 15) 

a 

Our lattice version of ^-functions are 

3 3 
6(x lt ) = a- 1 6 Xll0 , 5(x) = H S(x k ) , 6(x) = l[S(x,) (A.16) 

fc=l ^=0 

and we use 

6»(x M ) = 1 for Xfl > (A. 17) 

0{Xfj,) = otherwise 

Fields in momentum space are introduced by the Fourier transformation 

. \ f{x) = J h r Y. p e lpx Kp) in a T x i 3 volume 

I /O) = U /a (2Fpr e * pa 7(p) in infinite volume 
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A. 1.4 Continuum gauge fields 

An SU(A r ) gauge potential in the continuum theory is a vector field A^(x) with values 
in the Lie algebra su(N). It may thus be written as 

A M (x) = A«(x)T a (A.19) 

with real components A^(x) and 

(T a )t = -T°, tr {T a T b } = -\5 ab . (A.20) 
The associated field tensor, 
Ffjn/ {x 

) = dfj,A v (x) - d v An(x) + [Afj,(x), A v (x)\, (A. 21) 

may be decomposed similarly and the right and left action of the covariant derivative 
is defined by 

D^{x) = {d tl +A^{x), (A.22) 
^(x)D^=^(x)(d^~A^). (A.23) 

We note that periodic boundary conditions up to a phase 0^ are equivalent to adding 
a constant abelian gauge field iO^/L: in the above we replace A^ — > A^ + iO^/L. 

A. 1.5 Lattice action 

Let us first assume that the theory is defined on an infinite lattice. A gauge field U 
on the lattice is an assignment of a matrix U(x,ii) £ SU(iV) to every lattice point 
x and direction u = 0,1,2,3. Quark and anti-quark fields, ^p{x) and ip(x), reside on 
the lattice sites and carry Dirac, colour and flavour indices. The (unimproved) lattice 
action is of the form 

S[U^,^] = S G [U] + S F [U,^], (A.24) 

where Sq denotes the usual Wilson plaquette action and S F the Wilson quark action. 
Explicitly we have 

S vlU] = 4 £ tr {1 - U(p)} = 4 E E P ^ x ) > (A.25) 

P llv (x) = U(x, pi) U(x + a/t, v) U(x + ai>, y)^ 1 U(x, v)~ l (A.26) 

with go being the bare gauge coupling and U{p) the parallel transporter around the 
plaquette p. The sum runs over all oriented plaquettes p on the lattice, i.e. indepen- 
dently over fi, v. The quark action, 

S F [U,^,tP] =a 4 J2^(x)(Dw+mo)ip(x), (A.27) 

X 

is defined in terms of the Wilson-Dirac operator 

D w = \ { 7 „(V* + V M ) - aV* VJ , (A.28) 



which involves the gauge covariant lattice derivatives and V M , eq. (A. 9), and the 
bare quark mass matrix, m — diag(m 0u , mod, • • ■) ■ 
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A.1.6 Renormalization group functions and invariants 

Our RG functions are denned through 

1 d/i 
li dm 
fn d\i 

fj, a$ 



r(9) , 



(A.29) 
(A.30) 
(A.31) 



in terms of running coupling and running quark mass as well as some matrix ele- 
ment $ of a (multiplicatively renormalizable) composite field. They have asymptotic 
expansions 



9-S-O 



'f {bo + fh + . . .} 



bo 



(4tt) 



^(H-fJVf) , h 



r(9) -?{do+?d l + ...} 



do = 8/(4tt) 



102- 3v 

2 



7(5) 



9^0 



-5 2 {7o+.9 2 7i + ---} 



(A.32) 

(A.33) 
(A.34) 



The integration constants of the solutions to the RGEs define the RG invariants 



A = Ai (6 .g 2 )- bl/(2fc ° ) e- 1 /( 2 ^ 2 )exp 
M = fn (2b g 2 y do/2bo exp j- £ dx 
$ RGI = $[2fe .g 2 ]^ o/2fc0 



da; 



(x) ' b x 3 bix 



exp 



dx 



rjx) _ do^ 
f3(x) bog 

' l{x) 7o 



/3(a;) box 

where g = g(/j,) ... <f> = We will also us the shorthand notation 



A 

fj, 


= <Pg(S) 


= exp{- 


M 

fn 


= fm{g) 


= exp j- 


$RGI 

$ 


= 


= exp |- 



,(A.35) 
(A.36) 
(A.37) 

(A.38) 
(A.39) 
(A.40) 



dx^ 

with the constants exactly as defined above. 

A. 2 Conversion functions and anomalous dimensions 

Conversion functions and the anomalous dimensions 7 ma tch are not part of the stan- 
dard phenomenology literature. For completeness we give the explicit relations to the 
matching coefficients found directly in the literature and discuss the accuracy of their 
perturbative expansion. 
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A. 2.1 Matching coefficients and anomalous dimension 



We here describe the result (Bekavac et al, 2010) and its relation to the anomalous 



dimension. We denote a matrix element of some heavy-light quark bilinear ipTip^ in 
the effective theory by The Dirac structure T is left implicit Q 

All quantities are renormalized in the MS-scheme, with a scale \i a for the QCD 
bilinear and a scale /i in HQET. Choosing the pole quark mass toq|^] the matrix 
element is then (without explicit superscripts "QCD" we refer to HQET quantities, in 
the static approximation), 

$ QCD (m Q ,^ ;V kin ) = C ma t c h(m Q ,/i ,/i) x V kin ) + 0(l/m) . (A.41) 

The kinematical variables entering the matrix element $ are denoted by Vkin- For 
the (parti ally) conserved currents V^,A^ there is no /x G -dependence on the l.h.s. of 



eq. (A.41 1, (d fJlo <& cjcd (toq, fi ) = 0), while in general we have 

H 9$Q CD (m Q , Mo ) dln($Q CD (m Q , Mo )) _ 
^(m Q ,, o) Wo = dMVo) = %(3(Alo)) ' (A - 42) 

We pass to the RGI matrix element in QCD via (0(l/m) is dropped without notice) 

^gi D = exp {~ J 9(MO) dx J^} $ QCD (™Q.Mo;V kin ) (A.43) 

[ rSipa) lo ( x ) ] _ 
= expi - / dx— — } C matc h(TOQ,/i ,^) x Vkin) ■ (A.44) 
{ J J 

(0(l/m) is dropped without notice). It depends on the quark mass but not on a 
renormalization scale. The physical anomalous dimension is given by 

dln(m Q ) 91n(a match (m Q ,// 0) /i)) . 

7match(S*) = -TT-, \ ^ , (A.45) 

dm(m*J am(mQ) 



where the first factor is computed from the expansion ( 


Gray et al., 1990 


Fleischer 


et al., 1999 


Melnikov and Ritbergen, 2000 


Bekavac et al, 2010) 



-2 / \ 

m Q = m* [1 + X fcj[a(m*)]'] , , a{p) = (A.4 

fci = 4/3 , k 2 = -1.0414(JV t - 1) + 13.4434 , 
fc 3 = 0.6527(N f - l) 2 - 26.655(N f - 1) + 190.595 . 



1 The notation Cp of I Broad hurst and Grozin, 199 51) translates to our T a s T = (1, 70, 71, 7071) 



r = (75, 7075, 7fe, 7o7fc) ana ( |HeKavac et at., ^Uiu[ l uses the notation of proadhurst and Grozin^] 
|1995| when one sets d ( ,7 (1 = 70, 7± = Ik as it is the case in the rest frame. We will also refer to the 
bilinears as (PS, Aq, Vfc, T). In comparison to pekavac et al., 2010[) we add a subscript Q to the pole 
quark mass and a bar to the running mass (m — > niQ, m(/t) — ¥ for clarity. 

2 While in the complete, non-perturbative theory, the pole mass is ill-defined, in perturbation 
theory it exists order by order in the expansion. We use it here, because the formulae in the literature 
are written in terms of it. It will be eliminated in the final formulae. 
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The authors of Ref. (Bekavac et al, 2010) set fi a = fi. Building on ( Ji and Musolf, 1991 



Broadhurst and Grozin, 1995 Gimenez, 19921, they give the perturbative expansion 



C ma tch(TOQ,jU,M) = 1 + X! ^ L ik[M m Q/ l j2 )] k [«( m Q) 



Z>1 fe=0 



(A.47) 



with coefficients depending on the Dirac-structure, T. 
Independence of the l.h.s. of eq. ( |A.41 ) of fi yields 



dln(C match (TOQ, n ,n)) dln($(/x)) 
dln(n) <91n(/z) 



-7 stat (5(M)), 



(A.48) 



and with ain(5m gT n h ( ( "?^°' M)) = 7o(5(Mo)) we have 



d ln(C mat ch (wq , "t Q , m Q ) ) 
dln(mQ) 

91n((C ma tch(TOQ,/i ,/i)) 



<91n(m Q ) 



(A.49) 



+ lo{g(m Q )) - 7 sta t(ff(m Q )) . 



From these equations 7match(ff*) can be determined up to three- loop order and the 
differences TmatchO*) ~ 7match(s*) U P to four-loop order. 

A. 2. 2 Numerical results and the behavior of perturbation theory 

Let us now look at the numerical size of the perturbative coefficients of the RG func- 
tions. The following table lists results for Nf — 3. This is enough to understand the 
general picture since for smaller Nf the higher order coefficients are generically some- 
what larger, but not by much. 

coefficient i = 1 i = 2 i = 3 i = 4 



(4 7 r) i 6 i _ 1 
(4tt)H-i 
(47r) l 7 stat ,i-i 



0.71620 0.40529 0.32445 0.47367 
0.63662 0.76835 0.80114 0.90881 
-0.31831 -0.26613 -0.25917 



(^) i 7:i, i -i 

(4T) 4 7 I 7 natch, 4 -l 



~7fc 1 
'match, i— 1J 

'Y 75 ' 1 
"match, i— U 



-0.31831 
-0.31831 




-0.57010 
-0.87406 
0.30396 



-0.94645 
-3.12585 
2.17939 



14.803 



The normalization (An) 1 has been inserted such that the series is well behaved for 
a < 1/3 if the coefficients are order one. Indeed this is the magnitude of the coefficients 
in the first three rows which show as a comparison the beta-function, mass anomalous 
dimension and the anomalous dimension of the static-light bilinears (all in the MS- 
scheme). In contrast in the physical anomalous dimension of the vector current 7Z atc j 1 ) 
the 3-loop coefficient is rather big and the difference (47r) J [7match 3 — 7match 3] 1S even 
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above ten. Perturbation theory is then useful only at rather small a; in particular not 
really for the b-quark. 

An attempt to improve the perturbative series is to re-expand 7 ma tch in the coupling 
at a different scale, adjusting the scale to obtain smaller coefficients. In fact, since the 
effective theory is valid at energy scales below the mass of the quark, it is plausible that 
scales smaller than m* are more suitable. So we choose a coupling g 2 = g 2 (s^ 1 m i< ) = 
a{g 2 ,s) and 



7match(5) = 7matchO(3 2 , 1/«)] 1/2 ) , 

which is of course expanded order by order, 

g 2 = a(g 2 ,l/ S )=g 2 -2b ln(s)g 4 + .. 
The conversion functions are then expressed as 



C PS (M/A) = exp • 



rS , 7match(£) 

da;- 



(A.50) 



(A.51) 



(A.52) 



The difference comes from truncating eq. (A.50) as a series in g . The argument above 



suggests s > 1. The perturbative coefficients are listed in the following table for a few 
choices of s, for example the one which brings the two- loop coefficient 71 to zero. 



coefficient 


i = 1 


i = 2 


i = 3 


i = 4 


s 


(4^7^-1 


-0.31831 
-0.31831 


-0.57010 



-0.94645 
0.39720 




1 

3.4916 


(4^7^-1 


-0.31831 
-0.31831 


-0.87406 



-3.12585 
-0.231121 




1 

6.8007 


V ±n } L /match, i-l 'match, t-lJ 








0.30396 
0.30396 
0.30396 
0.30396 


2.17939 
0.972221 
-0.05414 
-0.23495 


14.803 
4.733 
1.82678 
1.85344 


1 

4 
13 
16 



The higher order coefficients can indeed be reduced significantly but s^4 is required. 
For B-physics airn^/s) is then not small and there is no really useful improvement 
for phenomenology, see Fig. |3.2| We emphasize, however, that with s w 4 the series is 
much better behaved for masses that are a factor two or more higher than the b-quark 



mass. The pattern visible in in the tables reflects itself in Fig. 3.2 

Let us finally mention that the same behavior is found for Cmatch(wQ, toq, toq) 
for all Dirac structures of the currents. Their perturbative expansion in a coupling 
g(mq/s) is better behaved for s^4 than for s = 1. 
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